NAVAL  POSTGRADUATE  SCHOOL 

Monterey ,  California 


THESIS 


FURTHER  STUDIES  IN  FILMWISE 

CONDENSATION  OF  STEAM  ON 

HORIZONTAL  FINNED  TUBES 

by 

Keith  Andrew  Swensen 

March,  1992 


Thesis  Co-Advisor 
Thesis  Co-Advisor 


P.J.  Marto 
S.B.  Memory 


Approved  for  public  release;  distribution  is  unlimited. 


T258717 


Unclassified 


CURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No  0704-0188 


i    REPORT  SECURITY  CLASSIFICATION 

Unclassified 


1b    RESTRICTIVE  MARKINGS 


SECURITY  CLASSIFICATION  AUTHORITY 


>   DECLASSIFICATION /DOWNGRADING  SCHEDULE 


3     DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
is  unlimited. 


PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


5    MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


NAME  OF  PERFORMING  ORGANIZATION 

Naval  Postgraduate  School 


6b    OFFICE  SYMBOL 
(If  applicable) 

34 


7a    NAME  OF  MONITORING  ORGANIZATION 

Naval  Postgraduate  School 


ADDRESS  (City,  State,  and  ZIP  Code) 

Monterey,  CA  93943-5000 


7b    ADDRESS  (City,  State,  and  ZIP  Code) 

Monterey,  CA  93943-5000 


NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


8b    OFFICE  SYMBOL 
(If  applicable) 


9    PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


ADDRESS  (City,  State,  and  ZIP  Code) 


10    SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 
NO 


TASK 
NO 


WORK  UNIT 
ACCESSION  NO 


title  (include  Security  Class,fication)  FURTHER   STUDIES    IN    FILMWISE    CONDENSATION    OF    STEAM  ON 

HORIZONTAL   FINNED   TUBES    (Unclassified) 


PERSONAL  AUTHOR(S) 


Keith   Andrew   Swensen 


a    TYPE  OF  REPORT 

Master's  Thesis 


13b    TIME  COVERED 
FROM  TO 


14.  DATE  OF  REPORT   (Year,  Month,  Day) 

March   1992 


15    PAGE  COUNT 
142 


supplementary  notation       The   views   expressed    in    this  thesis   are    those    of    the   author   and   do 
not    reflect    the    official    policy   or   position   of   the   Dept .    of   Defense    or   U.S.    Government 


cosati  codes 


FIELD 


GROUP 


SUB-GROUP 


18    SUBJECT  TERMS  {Continue  on  reverse  if  necessary  and  identify  by  block  number) 

filmwise  condensation,  integral  finned  tubes, 
vapor  velocity 


ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Over  the  years,  there  has  been  significant  variation  in  the  filmwise  steam 
condensation  data  at  NPS  on  horizontal  low- integral  finned  tubes.  With  a  view  to 
increasing  the  accuracy  of  the  data,  inserts  were  used  inside  the  tubes  to  reduce 
inside  thermal  resistance;  however,  significant  discrepancies  then  occurred  in  the 
calculated  outside  heat-transfer  coefficient  when  compared  to  data  taken  without  an 
insert.  These  discrepancies  arose  due  to  the  data  reduction  technique  which  assumes 
a  known  inside  heat-transfer  resistance  and  subtracts  this  from  a  measured  overall 
resistance.  If  the  assumed  value  on  the  inside  is  inaccurate,  then  the  outside  value 
is   equally  inaccurate. 

The  present  work  uses  an  instrumented  smboth  tube  to  obtain  accurate  inside  heat- 
transfer  correlations  both  with  and  and  without  inserts  and  uses  these  to  obtain 
accurate  outside  coefficients  for  a  family  of  uninstrumented  finned  tubes  with  a  view 
to  finding  an  optimum  fin  spacing  for  steam  condensation. 


DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
1)3  UNCLASSIFIED/UNLIMITED       D  SAME  AS  RPT  □  DTIC  USERS 


21    ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


a    NAME  OF  RESPONSIBLE  INDIVIDUAL 

P.J.    Marto 


22b  TELEPHONE  (Include  Area  Code) 

(408)    646-2989 


22c    OFFICE  SYMBOL 

69Mx 


'Form  1473,  JUN  86 


Previous  editions  are  obsolete 
S/N   0102-LF-014-6603 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

Unclassified 


Approved  for  public  release;  distribution  is  unlimited. 

Further  Studies  in  Filmwise 

Condensation  of  Steam  on 

Horizontal  Finned  Tubes 

by 

Keith  Andrew  Swensen 

Lieutenant,  United  States  Navy 

B.S.,  Brigham  Young  University,  1985 

Submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 

from  the 


ABSTRACT 

Over  the  years,  there  has  been  significant  variation  in  the  filmwise  steam 
condensation  data  at  NPS  on  horizontal  low-integral  finned  tubes.  With  a  view  to 
increasing  the  accuracy  of  the  data,  inserts  were  used  inside  the  tubes  to  reduce 
inside  thermal  resistance;  however,  significant  discrepancies  then  occurred  in  the 
calculated  outside  coefficient  when  compared  to  data  taken  without  an  insert.  These 
discrepancies  arose  due  to  the  data  reduction  technique  which  assumes  a  known 
inside  heat-transfer  resistance  and  subtracts  this  from  a  measured  overall  resistance. 
If  the  assumed  value  on  the  inside  is  inaccurate,  then  the  outside  value  is  equally 
inaccurate. 

The  present  work  uses  an  instrumented  smooth  tube  to  obtain  accurate  inside 
heat-transfer  correlations  both  with  and  without  inserts  and  uses  these  to  obtain 
accurate  outside  coefficients  for  a  family  of  uninstrumented  finned  tubes  with  a  view 
to  finding  an  optimum  fin  spacing  for  steam  condensation. 
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I.   INTRODUCTION 

A.       BACKGROUND 

A  reduction  in  size  and  weight  of  all  types  of  heat  exchangers  aboard  Naval 
vessels  will  allow  more  efficient  use  of  space.  The  benefits  might  include  greater 
equipment  accessibility  for  maintenance  or  greater  heat  exchanger  capacity  (without 
a  corresponding  increase  in  size  and  weight)  with  a  corresponding  increase  in  fuel 
efficiency. 

For  the  past  ten  years,  the  Naval  Postgraduate  School  in  collaboration  with  the 
David  Taylor  Research  Center  and  the  National  Science  Foundation  has  conducted 
research  that  is  directed  at  the  development  of  smaller,  more  efficient  steam 
condensers.  Improved  designs  can  result  in  significant  space  savings,  always  a 
primary  concern  on  Naval  vessels,  especially  submarines. 

Uncertainties  in  past  data  using  steam  were  apparently  due  to  the  lack  of 
detailed  information  about  the  inside  heat-transfer  correlations  used  during  the  data 
reduction  process.  Previously,  the  standard  Sieder-Tate  correlation  was  assumed  to 
be  valid  for  the  inside  heat-transfer  coefficient,  but  it  may  not  be  the  best  correlation 
to  use  with  the  particular  test  arrangement  used  in  this  research  program. 

A  large  amount  of  enhanced  condensation  data  has  been  collected  in  previous 
studies  at  the  Naval  Postgraduate  School  on  more  than  90  different  condenser  tubes 
of  varying  fin  height,  fin  spacing,  and  tube  material,  with  a  view  to  finding  an 


optimum  fin  geometry  for  both  steam  and  refrigerant  condensation.  Satisfactory 
results  have  been  obtained  with  refrigerant  data  using  R-113.  However,  some 
troublesome  questions  of  possible  contamination  and  instrument  inaccuracy  still 
remain  with  the  steam  data.  It  is  felt  that  in  order  to  address  these  questions,  a 
fundamental  evaluation  of  the  heat  transfer  apparatus  on  which  this  data  was 
collected  and  of  the  data  reduction  process  was  considered  appropriate.  Current 
data,  taken  on  a  carefully  cleaned  and  calibrated  apparatus,  could  be  compared  to 
previously  recorded  data  and  a  determination  as  to  its  validity  and  reproducibility 
could  be  made.  In  addition,  a  thorough  evaluation  of  the  best  inside  heat-transfer 
correlation  would  lead  to  more  reliable  steam  condensation  results. 

B.       CONDENSATION 

Condensation  occurs  when  a  vapor  is  cooled  below  its  saturation  temperature, 
or  when  a  vapor/gas  mixture  is  cooled  below  its  dew  point.  Surface  condensation 
occurs  in  condensers  when  a  cooled  surface  (kept  at  a  temperature  below  the 
saturation  temperature  of  the  vapor)  contacts  the  vapor.  The  vapor  molecules  that 
contact  such  a  surface  stick  to  that  surface  and  condense  into  liquid  molecules. 
Condensation  may  occur  in  one  of  the  following  modes:  filmwise,  dropwise,  or  mixed 
mode  (a  combination  of  filmwise  and  dropwise)  condensation.  In  the  filmwise  mode, 
the  liquid  wets  the  cold  surface  to  form  a  continuous  film.  If  the  liquid  does  not  wet 
the  surface  but  instead  forms  discrete  drops  on  the  cold  surface,  dropwise  or  mixed 


mode  condensation  will  occur  and  is  often  caused  by  some  form  of  contamination. 
[Ref.  1] 

The  condensate  forming  on  the  tube  surface  offers  a  resistance  to  heat  transfer 
between  the  vapor  and  the  surface,  which  increases  with  the  thickness  of  the  liquid 
layer.  Even  though  dropwise  condensation  results  in  much  larger  heat-transfer 
coefficients  than  filmwise  condensation,  it  is  difficult  to  maintain  a  stable  dropwise 
condition  over  prolonged  periods.  Therefore,  in  most  cases  condenser  design 
calculations  are  based  on  the  assumption  of  filmwise  condensation,  resulting  in  lower 
heat-transfer  coefficients  and  more  conservative  designs.  [Ref.  1] 

In  a  condenser,  the  coolant  side,  tube  wall,  and  vapor  side  thermal  resistances 
control  the  heat  transfer  rate  from  vapor  to  coolant.  Also,  for  experimental  work 
we  always  use  clean  tubes,  but  in  real  condensers  tubeside  fouling  can  play  an 
important  role  an  increasing  the  coolant  side  resistance.  The  magnitudes  of  these 
resistances  depend  on  the  fluid,  tube  geometry,  and  flow  conditions  on  the  vapor  and 
coolant  side.  For  steam  condensation,  it  is  the  coolant  side  thermal  resistance  which 
tends  to  dominate.  Methods  to  lower  this  inside  resistance  include  the  use  of  inserts 
or  roped  tubing  to  promote  turbulence,  thereby  raising  the  convection  heat-transfer 
coefficient.  However,  such  modifications  lead  to  increased  pressure  drop  through 
the  tubes,  which  must  be  compensated  for  by  providing  extra  pumping  capacity. 
Heat  transfer  through  the  tube  wall  is  conductive  and  is  fixed  once  tube  thickness 
and  material  are  selected.  The  vapor  side  resistance  is  due  to  the  condensate  film 
which  forms  on  the  outside  of  the  tube.    For  filmwise  condensation,  the  outside 


resistance  can  be  reduced  by  the  addition  of  low  integral  fins.  These  have  the  effect 
of  not  only  increasing  the  outside  surface  area  of  the  tube,  but  also  of  thinning  the 
condensate  film  around  the  fins  due  to  surface  tension  forces.  Too  small  a  fin 
spacing  may  result  in  condensate  flooding,  whereas  too  large  a  fin  spacing 
approaches  the  smooth  tube  case;  there  should  be  an  optimum  fin  spacing 
somewhere  in  between  these  two  extremes.  Horizontal  fin  spacing  is  therefore  of 
prime  importance,  and  finding  the  optimum  spacing  is  one  of  the  objectives  of  this 
long-term  research  program. 

C.      CONDENSATION  RESEARCH  AT  NAVAL  POSTGRADUATE  SCHOOL 

The  research  effort  at  NPS  has  included  the  study  of  differing  fin  dimensions 
(i.e.  fin  height,  fin  width,  fin  spacing)  on  low-integral  finned  horizontal  tubes. 
Experimentation  has  included  the  use  of  three  different  test  fluids  (steam,  R-113, 
and  ethylene  glycol)  under  various  operational  conditions  using  a  number  of  different 
tube  diameters. 

Van  Petten  [Ref.  2]  provides  a  summary  of  research  efforts  at  NPS  through  the 
end  of  1988.  Van  Petten  and  subsequent  researchers  have  analyzed  small,  medium, 
and  large  diameter  finned  tubes  to  find  the  optimum  fin  spacing  for  maximum  heat- 
transfer  enhancement  of  the  fluids  mentioned  above.  However,  discrepancies  found 
by  Guttendorf  [Ref.  3]  in  the  data  processing  technique  (modified  Wilson  plot), 
which  resulted  in  different  values  of  heat-transfer  enhancement  (for  the  same  tube 


under  the  same  operating  conditions  depending  on  whether  an  insert  was  or  was  not 
used),  have  raised  doubts  about  the  accuracy  of  the  inside  heat- transfer  correlation. 
Rouk  [Ref.  4]  investigated  the  use  of  an  optimization  technique  to  predict  the 
inside  heat-transfer  correlation.  When  the  optimization  effort  proved  unsuccessful, 
he  next  used  the  instrumented  smooth  tube  data  of  Georgiadis  [Ref.  5]  to  develop 
an  inside  heat  transfer  correlation,  but  could  not  find  a  correlation  with  sufficient 
accuracy  based  on  previous  data.  He  recommended  that  once  an  overhaul  on  the 
test  apparatus  was  complete,  an  increase  in  data  precision  would  allow  the 
development  of  an  accurate  inside  heat-transfer  correlation.  This  work  is  a  follow 
on  effort  to  develop  inside  heat-transfer  correlations  which  can  predict  the  value  of 
the  inside  heat-transfer  coefficient  with  good  accuracy  under  a  variety  of  flow 
conditions.  Once  an  inside  correlation  can  be  found,  the  object  of  this  effort  is  to 
reprocess  previous  data  and  see  if  the  discrepancies  reported  by  former  researchers 
for  finned  tubes  on  this  apparatus  can  be  rectified. 

D.      OBJECTIVES 

The  main  objectives  of  this  thesis  were  to: 


1.  Disassemble  and  meticulously  clean  the  apparatus  to  eliminate  any  existing 
contamination  with  a  view  to  eliminating  dropwise  condensation  problems 
experienced  in  the  past. 

2.  Carefully  reassemble  the  apparatus  using  new  gasket  material,  and  make 
modifications  to  improve  system  performance. 

3.  Recalibrate  all  system  instrumentation  to  ensure  the  greatest  achievable 
accuracy. 


4.  Investigate  the  possibility  of  manufacturing  large,  medium,  and  small  diameter 
instrumented  smooth  tubes. 

5.  Use  the  new  instrumented  tubes  and  the  one  existing  medium  diameter 
smooth  instrumented  tube  (fabricated  by  Poole  [Ref.  6])  to  obtain  accurate 
inside  heat-transfer  correlations  for  a  number  of  insert  types  as  well  as  the 
no  insert  condition. 

6.  Evaluate  the  accuracy  of  the  currently  used  data  processing  technique 
(modified  Wilson  plot)  using  instrumented  tube  data. 

7.  Reprocess  previous  data  using  the  new  correlations  with  a  view  to  comparing 
current  and  past  smooth  tube  and  finned  tube  data  to  provide  continuity  with 
previous  studies. 


II.    LITERATURE  SURVEY 

A.       INTRODUCTION 

When  a  vapor  condenses  in  the  filmwise  mode  on  a  smooth  horizontal  tube  it 
forms  a  thin  continuous  film  of  condensate  on  the  surface  of  the  tube.  The 
condensate  film  thickens  around  the  tube  due  to  gravity.  This  condensate  film 
provides  a  resistance  to  heat  transfer  which  may  be  lowered  through  the  use  of  fins. 
For  quite  some  time,  it  was  thought  to  be  impractical  to  use  finned  tubes  with  high 
surface  tension  fluids  such  as  water,  due  to  condensate  retention  and  flooding 
between  the  fins.  However,  a  number  of  studies  conducted  on  finned  tubes  using 
steam  have  shown  that  substantial  heat-transfer  enhancement  may  be  achieved. 

A  significant  amount  of  research  at  the  Naval  Postgraduate  School  and 
elsewhere  has  addressed  the  issue  of  optimum  fin  height,  thickness,  and  spacing 
required  for  maximum  heat  transfer.  Yau  et  al  [Ref.  7]  reported  that  "with  an 
increase  in  fin  density,  up  to  a  limit  (this  limit  is  not  yet  known  in  a  generalized 
manner),  the  heat-transfer  coefficient  increases  at  a  rate  faster  than  the  increase  in 
the  outside  area  due  to  the  presence  of  fins.  This  additional  enhancement  is  due  to 
the  thinning  effect  of  the  surface-tension  forces  on  the  condensate  film. 
Unfortunately,  surface-tension  forces  also  adversely  affect  heat  transfer  by  causing 
condensate  to  be  retained  between  fins"  [Ref.  8].  Katz  et  al  [Ref.  9]  also  found  that 


on  finned  tubes  the  portion  of  the  surface  occupied  by  condensate  is  dependent  upon 
the  ratio  of  condensate  surface-tension  to  density  and  the  fin  geometry. 

Condensate  retention  and  the  behavior  of  the  condensate  film  on  the  tube 
surface  under  various  conditions  are  critical  parameters  in  the  heat  transfer  process 
on  horizontal  finned  tubes.  Several  models  have  been  developed  to  predict  this 
behavior  and  the  reader  is  referred  to  an  extensive  review  of  horizontal  finned  tube 
heat  transfer  by  Marto  [Ref.  10]  for  a  more  detailed  coverage  of  the  topic. 

B.      VAPOR  SIDE  CONSIDERATIONS 

The  filmwise  condensation  of  vapor  on  a  horizontal  tube  is  a  complex  two- 
phase  heat  transfer  process,  for  which  a  suitably  complex  model  would  be  required 
to  accurately  predict  heat  transfer  performance  under  all  conditions. 

In  1916,  Nusselt  [Ref.  11]  set  forth  his  theoretical  work  on  the  study  of  laminar 
filmwise  condensation  of  a  "stationary"  vapor  on  a  vertical  or  inclined  plate  and  a 
horizontal  tube.   Nusselt's  simplifying  assumptions  included  the  following  [Ref.  12]: 

1.  Pure  saturated  vapor 

2.  Negligible  vapor  velocity  (11^=0) 

3.  Heat  transfer  across  the  condensate  film  by  conduction  only 

4.  Laminar  condensate  flow  governed  only  by  gravitational  and  viscous  forces 

5.  Condensate  properties  constant 

6.  Isothermal  condensing  surface 

7.  Negligible  interface  temperature  drop 

8 


Nusselt's  result  for  the  mean  heat-transfer  coefficient  for  a  horizontal  tube  was 


obtained: 


Nu  =0.728 


*/P/P/-Pjg*» 


H/D.Ar, 


1/4 


(2.1) 


or 


where: 


Nu 

kf 

Pf 

Pv 

g 

\ 

AT, 


Nu  =0.655 


*/P/P/-Pv)g^ 


1/3 


(2.2) 


mean  Nusselt  number 

thermal  conductivity  of  condensate  film  (W/m  •  k) 

condensate  film  density  (kg/m3) 

vapor  density  (kg/m3) 

gravitational  constant  (9.81  m/s2) 

specific  enthalpy  of  vaporization  (J/kg) 

dynamic  viscosity  of  condensate  film  (N-s/m2) 

outside  tube  diameter  (m) 

average  temperature  difference  across  condensate  film  (K) 

heat  flux  based  on  outside  area  (Q/A0)  (W/m2) 


Many  workers  have  improved  on  Nusselt's  theoretical  analysis  by  accounting 
for  some  of  the  terms  he  neglected  through  his  simplifying  assumptions.  However, 
equations  (2.1)  and  (2.2)  have  been  found  to  be  remarkably  accurate  over  a  wide 
range  of  conditions  for  a  stationary  vapor.  High  vapor  velocity  can  increase  film 
condensation  heat  transfer  substantially.  This  enhancement,  which  refers  to  the 
amount  of  heat  transfer  above  or  below  the  value  predicted  by  the  Nusselt  analysis, 
is  due  to  the  effect  of  thinning  the  condensate  film.  However,  vapor  shear  is  the  one 


assumption  which  if  applied  can  lead  to  significant  increases  in  the  heat-transfer 
coefficient. 

The  theoretical  result  of  Shekriladze  and  Gomelauri  (1966)  [Ref.  13],  who 
considered  interfacial  shear  stress  due  to  vapor  velocity,  is  shown  in  equation  (2.3). 

=0.64(1  +(1  +  1.69f)1/2)1/2  (2.3) 


Ke2* 


where: 


Nu     =    Nusselt  number  for  the  vapor  side 

Re24>  =    two  phase  Reynolds  number,  (pf  U.  DJfif) 


For  steam  condensation,  the  empirically  derived  correlation  of  Fujii  et  al  [Ref. 
14]  is  shown  in  equation  (2.4). 

=  0.96F^  (2.4) 


Re2<St 


The  Nusselt  expression  (equation  (2.1)  can  be  expressed  in  similar  form: 


^L=0.728F1/4  (2.5) 


Re2* 


Whereas  the  vapor  velocity,  U.,  cancels  out  in  the  Nusselt  expression 
(stationary  vapor  assumption),  the  Fujii  correlation  includes  the  vapor  velocity  effect. 
Therefore  we  can  expect  equation  (2.4)  to  more  accurately  predict  steam-side  heat 
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transfer  coefficients  for  those  cases  where  vapor  velocity  begins  to  have  a  significant 
impact. 

For  further  review  of  basic  theoretical  studies  on  the  subject  of  laminar  film 
condensation  on  smooth  tubes  the  reader  is  referred  to  Rose  [Ref.  15]. 

C.       COOLANT  SIDE  CONSIDERATIONS 

For  a  turbulent  flow  regime  inside  a  pipe  (Re  >  10,000),  a  number  of  coolant- 
side  correlations  have  been  used;  many  of  these  have  taken  the  form: 

Nu=CiRemPrn  (2.6) 


where: 


Nu  =  mean  coolant  Nusselt  number  for  turbulent  flow 

C,  =  correlating  coefficient 

Re  =  coolant  Reynolds  number 

Pr  =  coolant  Prandtl  number 


The  most  common  correlations  with  the  same  form  as  equation  (2.6)  are  that 
of  Dittus  and  Boelter  (1930)  [Ref.  16]: 


Nu  =0.023  Re°*Pr0A  (2J) 


and  Colburn  (1933)  [Ref.  17]: 

Nu  =0.023  Re^Pr113  (2-8) 
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Sieder  and  Tate  (1936)  [Ref.  18]  applied  a  correction  factor  to  equation  (2.8) 
to  account  for  cases  in  which  the  bulk  to  inner  wall  temperature  difference  is  large 
enough  to  cause  substantial  variations  in  coolant  viscosity  as  follows: 

Nu  =  0.027  Re0*Prl/*(\iJ\i  J  u  (2.10) 

where: 

/xc    =    coolant  viscosity  evaluated  at  mean  bulk  temperature  (N-s/m2) 
/iw    =    coolant  viscosity  evaluated  at  mean  inner  tube  wall  temperature 

(N-s/m2) 

The  fluid  properties  in  equations  (2.7),  (2.8),  and  (2.9)  are  evaluated  at  mean 
coolant  bulk  temperature  Tmbulk. 

T*u,-V,*TJP  (2.11) 

where: 

Tj       =    tube  coolant  inlet  temperature  (K) 
T2      =    tube  coolant  outlet  temperature  (K) 

Equations  (2.7),  (2.8)  and  (2.9)  are  valid  for  Re>  104  and  0.7  <  Pr  <  100,  and 
were  developed  for  long  smooth  pipes  with  no  inserts  [Ref.  12]. 

The  use  of  inserts  and  the  effect  of  bends  close  to  the  tube  entrance  region 
can  affect  the  values  of  both  the  leading  coefficient,  C,  ,  and  Reynolds  number 
exponent,  m,  in  equation  (2.6).     One  of  the  major  focuses  of  this  study  is  to 
determine  the  values  of  C,  and  m  for  the  no  insert,  wire  wrap  insert,  and  Heatex 
insert  cases. 
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Other  well-known  turbulent  pipe  flow  heat-transfer  correlations  (i.e.  Petukhov- 
Popov  [Ref.  19],  Sleicher-Rouse  [Ref.  20],  etc.)  and  the  results  of  an  ANL  (Argonne 
National  Laboratory)  [Ref.  21]  study  which  evaluated  several  such  correlations  for 
accuracy  are  reviewed  in  section  VI  C. 
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III.  APPARATUS  AND  SYSTEM  INSTRUMENTATION 

A.       SYSTEM  OVERVIEW 

The  apparatus  used  for  this  research  was  basically  the  same  as  was  used  by  Van 
Petten  [Ref.  2]  and  Guttendorf  [Ref.  3]  with  certain  modifications.  A  system 
schematic  is  provided  in  Figure  1.  Steam  generated  from  distilled  water  in  the  .30 
m  diameter  pyrex  glass  boiler  using  ten  4  kW,  440  V  Watlow  immersion  heaters  was 
the  working  medium  for  this  set  of  experiments.  From  the  boiler  section  the  steam 
passed  up  through  a  reducing  section  and  a  2.13  m  straight  length  of  pyrex  glass 
piping,  (ID  of  0.15  m),  it  was  then  turned  through  180  degrees  using  two  90  degree 
pyrex  glass  elbows,  and  then  descended  down  a  1.52  m  straight  length  of  pyrex  glass 
piping.  The  steam  then  entered  the  stainless  steel  test  section  containing  the 
horizontally  mounted  condenser  tube  (see  Figures  1  and  2);  any  steam  not 
condensing  there  was  condensed  in  the  auxiliary  condenser  located  just  beneath  the 
test  section.  The  auxiliary  condenser  was  constructed  of  a  single  copper  coil  mounted 
to  a  stainless  steel  baseplate  enclosed  within  a  pyrex  glass  condenser  section. 
Coolant  flow  through  the  auxiliary  condenser  was  used  to  control  system  pressure, 
and  all  condensed  liquid  was  returned  via  the  condenser  baseplate  drain  to  the  boiler 
section  by  gravity. 

Coolant  for  the  auxiliary  condenser  was  provided  via  a  throttled  water 
connection  with  associated  flowmeter.    Coolant  flow  through  the  single  horizontal 
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Figure  1.   Schematic  of  Single  Tube  Test  Apparatus 
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Figure  2.   Schematic  of  Test  Section  Insert 
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tube  was  provided  by  a  separate  system  consisting  of  a  sump  tank  with  two 
centrifugal  pumps  connected  in  series.  Coolant  flow  rate  was  measured  by  a 
carefully  calibrated  flowmeter.  By  varying  coolant  flow-rate  through  the  single 
horizontal  tube,  the  rate  of  steam  condensation  on  the  tube  (and  hence  heat-transfer 
coefficient)  could  be  varied. 

Non-condensible  gases  were  removed  using  the  vacuum  pump  system  shown  in 
Figure  3.  The  condensing  coil  for  this  purge  system,  located  in  the  sump  tank, 
served  to  condense  steam  carried  through  the  vacuum  line  during  the  purging 
process.  The  vacuum  line  took  its  suction  from  the  base  of  the  auxiliary  condenser, 
the  coolest  spot  in  the  apparatus  and  the  place  where  non-condensible  gases  (i.e.  air) 
were  most  likely  to  accumulate. 

B.       SYSTEM  INSTRUMENTATION 

The  power  to  the  440  V  heaters  was  controlled  through  a  panel  mounted 
potentiometer.  A  description  of  the  power  calculation  for  input  into  the  data 
acquisition  system  can  be  found  in  Poole  [Ref.  6]. 

System  pressure  was  monitored  in  three  ways: 

1.  A  Setra  model  204  pressure  transducer 

2.  A  Heise  solid  front  pressure  gauge  (visual  reading  only) 

3.  System  saturation  pressure  from  vapor  temperature  measurement 
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Figure  3.   Schematic  of  Purging  System  and  Cooling  Water  Sump 


18 


System  vapor  temperature  was  monitored  using  both  teflon,  and  metal  sheathed 
type-T  copper/constantan  thermocouples  located  juxtaposed  in  the  test  section;  this 
position  was  just  upstream  of  the  test  condenser  tube.  Condensate  temperature  was 
also  monitored  using  a  teflon  coated  type-T  copper/constantan  thermocouple  located 
on  the  condensate  return  line  between  the  auxiliary  condenser  and  boiler.  Coolant 
temperature  rise  in  the  condenser  tube  was  measured  using  four  methods: 

1.  Two  teflon  coated  type-T  copper/constantan  thermocouples 

2.  Two  metal  sheathed  type-T  copper/constantan  thermocouples 

3.  Two  Hewlett-Packard  2804A  quartz  crystal  thermometers 

4.  A  ten-junction  teflon  coated  type-T  copper/constantan  thermopile 

These  were  all  placed  at  the  inlet  to  and  exit  from  the  condenser  tube;  at  the 
outlet,  all  thermocouples  were  placed  just  downstream  of  a  coolant  mixing  chamber. 

Two  data  reduction  programs  were  used  to  collect  and  reduce  data  on  this 
apparatus;  "DRPINST',  and  "DRPKS".  The  instrumented  tube  constructed  by  Poole 
[Ref.  6]  was  used  to  determine  an  accurate  inside  heat-transfer  correlation  for  inserts 
used;  this  instrumented  tube  contained  six  wall  thermocouples.  For  the  instrumented 
tube  the  appropriate  calibration  equations  were  accessed  in  the  data  acquisition 
program  "DRPINST".  For  non-instrumented  tubes,  the  data  reduction  program 
"DRPKS"  was  used.  Fluid  property  equations  used  in  the  data  reduction  programs 
are    given    in    Appendix   A    and    calibrations   were    conducted    for    all    system 
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instrumentation  (flowmeter,  thermocouples,  pressure  transducer,  etc.)  and  are 
included  in  Appendix  B. 

The  data  as  monitored  by  the  aforementioned  system  instrumentation  was 
processed  by  an  HP-3497A  data  acquisition  system  controlled  by  an  HP-9826A 
computer  provided  with  the  correct  data  acquisition  program.  The  raw  data  was 
processed  and  stored  on  computer  disks.  Program  channel  assignments  are  given  in 
Table  1. 

C.      SYSTEM  MODIFICATIONS 

At  the  beginning  of  this  investigation  the  apparatus  was  entirely  disassembled 
to  facilitate  complete  overhaul  of  the  system.  Modifications  and  details  of  assembly 
were  as  follows: 


1.  The  apparatus  was  taken  apart  piece  by  piece,  inspected,  meticulously  cleaned 
with  a  warm  solution  of  Sparkleen  biodegradable  soap  and  subjected  to  a 
complete  acetone  rinse  prior  to  reassembly. 

2.  A  new  pyrex  glass  riser  section  above  the  boiler  was  0.31  m  shorter  than  the 
previous  section  and  allowed  the  addition  of  a  new  aluminum  stand  on  which 
to  place  the  heater  baseplate.  This  new  stand  allowed  much  easier  access  to 
the  440  V  heater  wiring  plus  the  adjustable  legs  allowed  level  adjustment  of 
the  entire  apparatus,  to  ensure  proper  alignment  of  the  single  horizontal  tube. 

3.  The  two  pyrex  glass  elbows  were  replaced. 

4.  Every  gasket  in  the  system  was  replaced  (using  Buna-N  rubber)  and,  using  a 
standard  star  torque  pattern,  all  flanged  joints  were  tightened  to  a  final 
torqued  of  30  inch-pounds  (manufacturer  recommended  maximum  torque 
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Table  1.   DATA  ACQUISITION  SYSTEM  CHANNEL  ASSIGNMENT 

CHANNEL 

QUANTITY 
MEASURED 

APPLICABLE  TYPE 
DESIGNATION 

40 

Vapor  Temperature 

T-57 
(See  Appendix  B.l) 

41 

Vapor  Temperature 

T-56 
(See  Appendix  B.l) 

42 

Room  Temperature 

T-58 
(See  Appendix  B.l) 

43 

Tube  Coolant  in  1 

T-58 

44 

Tube  Coolant  out  1 

T-58 

45 

Tube  Coolant  in  2 

T-55 
(See  Appendix  B.l) 

46 

Tube  Coolant  out  2 

T-55 

47 

Condensate  Return 

T-58 

48 

Instrumented  Tube 

T-57  or  T-55 
as  applicable 

49 

Instrumented  Tube 

T-57  or  T-55 
as  applicable 

50 

Instrumented  Tube 

T-57  or  T-55 
as  applicable 

51 

Instrumented  Tube 

T-57  or  T-55 
as  applicable 

52 

Instrumented  Tube 

T-57  or  T-55 
as  applicable 

53 

Instrumented  Tube 

T-57  or  T-55 
as  applicable 

60 

10  -  Junction 
Thermopile 

T-61 
(See  Appendix  B.l) 

61 

Voltage 

NA 

62 

Current 

NA 

64 

Pressure  Transducer 

(See  Appendix  B.4) 
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specification  was  60  inch-pounds).   The  bolts  holding  the  flanges  together 
should  be  checked  for  tightness  on  a  regular  basis,  since  the  thermal  cycling 
of  the  apparatus  has  been  shown  to  cause  loosening  of  nut  and  bolt 
assemblies. 

5.  The  previous  vacuum  pump  was  replaced  with  a  Gast  model  2567-V108 
vacuum  pump  which  could  draw  vacuum  to  130  mmHg  with  an  installed  check 
valve  to  prevent  pump  back  spin  when  the  pump  was  stopped.  The  new 
pump  could  draw  vacuum  much  more  rapidly,  but  could  not  operate  at  less 
than  130  mmHg.  Remaining  non-condensible  gases  were  removed  by  flushing 
the  system  with  steam.  The  steam  flushing  procedure  for  removal  of  non- 
condensible  gases  is  given  in  the  operating  procedures  section.  Once  the 
system  is  completely  filled  with  steam,  operating  pressures  well  below  130 
mmHg  could  be  achieved  utilizing  the  auxiliary  condenser. 

6.  The  double  coil  auxiliary  condenser  was  replaced  with  a  single  coil  used 
originally  by  Van  Petten  [Ref.  2].  The  single  coil  was  not  coated  with  the 
special  oxide  coating  used  by  Guttendorf  [Ref.  3].  It  was  felt  to  be  superior 
to  the  double  coil  in  that  the  baseplate  welds  were  of  much  higher  quality  and 
were  preferred  on  the  basis  of  vacuum  tightness. 

7.  The  aluminum  side  plates  attached  to  the  pyrex  glass  auxiliary  condenser 
housing  were  replaced  with  new  stainless  steel  side  plates  with  penetrations 
for  pressure  bleed,  vacuum  line,  and  a  pressure  transducer.  These  three 
penetrations  were  fitted  with  screw  threaded  stainless  steel  connectors.  The 
stainless  steel  connectors  were  heli-arc  welded  in  place.  Prior  to  this 
modification  (completed  24  January  1992)  a  leak  test  conducted  from  21 
December  1991  to  2  January  1992  revealed  a  mean  vacuum  leak  rate  of  ~  3.4 
mmHg  per  day  (see  Figure  C.l).  A  subsequent  leak  test  conducted  from  6 
February  1992  to  19  February  1992  showed  an  improvement  in  the  mean  leak 
rate  to  —1.7  mmHg  per  day  (see  Figure  C.2). 

8.  System  instrument  modifications  included  the  addition  of  the  Setra  pressure 
transducer  and  the  Heise  pressure  gauge,  and  the  removal  of  the  mercury 
manometer.  All  system  instrumentation  was  recalibrated  and  the  results 
incorporated  into  the  data  reduction  programs. 

9.  Finally  the  apparatus  was  lagged  with  Halstead  insulating  foam  to  reduce  heat 
loss  as  much  as  possible.  The  test  section,  which  was  left  uncovered 
previously,  was  also  lagged. 
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IV  EQUIPMENT  OPERATION  AND  EXPERIMENTAL  PROCEDURE 

A.       SYSTEM  STARTUP  AND  SHUTDOWN  PROCEDURES 

Startup  of  the  system  is  accomplished  in  the  following  manner: 

1.  Ensure  distilled  water  level  in  the  boiler  is  4  to  6  inches  above  the  top  of  the 
heating  elements.  The  boiler  is  filled  by  gravity  drain  via  a  hose  connection 
from  the  distilled  water  tank  to  the  boiler  fill  valve.  Ensure  the  vent  valve  on 
the  side  of  the  auxiliary  condenser  is  open  when  filling  or  draining  the  boiler. 
The  boiler  may  be  drained  by  removing  the  hose  connection  and  opening  the 
fill  valve,  which  allows  drainage  into  the  trench  directly  beneath  the  boiler. 

2.  Once  the  boiler  is  filled  to  the  appropriate  level,  shut  the  boiler  fill  valve  and 
the  distilled  water  tank  valve. 

3.  Shut  the  system  vent  valve. 

4.  Turn  on  the  data  acquisition  system,  computer  and  printer.  Load  the 
appropriate  program  (DRPINST  or  DRPKS)  and  check  for  proper  operation. 
Then  check  all  thermocouple  outputs,  by  stepping  through  the  appropriate 
data  acquisition  system  channels,  to  verify  that  all  are  registering  ambient 
temperature. 

5.  Open  the  fill  valve  to  the  coolant  water  sump  tank  to  a  level  such  that  the 
tank  overflow  drain  box  does  not  overflow  (the  valve  is  located  between  the 
boiler  control  panel  and  heat  pipe  apparatus). 

6.  Turn  on  the  cooling  water  supply  pumps  and  adjust  the  tube  flow  rate  from 
20%  to  60%  of  the  rotameter  setting  and  check  for  leaks.  Reset  flow  rate  to 
desired  level. 

7.  Open  valves  from  tap  water  system  to  auxiliary  condenser  and  adjust  coolant 
flow  rate  to  at  least  30%  and  check  for  leaks.  Reset  flow  rate  to  at  least 
10%. 

8.  Energize  heaters  and  adjust  voltage  to  approximately  50  volts  (40  volts  if  the 
system  is  already  at  vacuum  below  100  mmHg  to  limit  the  vibrational  shock 
to  the  system  from  oversized  vapor  bubble  formation).  To  energize  the 
heaters  there  are  three  switches  which  must  be  placed  in  the  on  position.  The 
first  is  located  in  power  panel  p5  located  in  the  main  hallway  adjacent  to  the 
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lab  and  is  labeled  switch  3  /  heater  controller  room  106.  The  second  is  the 
heater  load  bank  circuit  breaker  located  on  the  side  of  the  boiler  control 
panel.  The  third  is  the  condensing  rig  boiler  power  switch  on  the  front  of  the 
boiler  control  panel.  Increase  the  voltage  gradually  in  10  volt  increments  to 
the  desired  level. 

9.  Turn  on  the  vacuum  pump  and  open  the  vacuum  line  valve.  Allow  the 
vacuum  pump  to  run  until  system  pressure  is  below  3  psi,  then  shut  the 
vacuum  line  valve  just  prior  to  turning  off  the  vacuum  pump. 

10.  As  system  warmup  continues  and  pressure  increases  to  above  4  psi  energize 
the  vacuum  pump  as  necessary  to  flush  the  non-condensible  gases  out  of  the 
system  through  the  vacuum  line  by  forcing  the  gases  out  with  steam.  To 
ensure  that  the  non-condensibles  gather  at  the  base  of  the  auxiliary 
condenser,  where  the  vacuum  line  suction  is  located,  ensure  that  the 
horizontal  tube  is  not  supplied  with  coolant  flow,  and  adjust  coolant  flow 
through  the  auxiliary  condenser  as  necessary  to  ensure  steam  is  filling  the 
entire  system.  The  auxiliary  condenser  may  be  touched  lightly  by  hand,  along 
its  entire  length,  to  ensure  the  system  is  completely  filled  with  steam;  any  cool 
spots  indicate  the  presence  of  non-condensible  gases  which  means  that  the 
flushing  process  is  not  complete.  The  flushing  process  takes  15  to  30  minutes 
to  accomplish,  and  should  be  repeated  periodically  for  long  periods  of 
operation. 

11.  At  the  conclusion  of  the  flushing  process,  shut  the  vacuum  line  valve  and 
secure  the  vacuum  pump. 

12.  In  order  to  ensure  that  filmwise  condensation  occurs  on  the  tube,  coolant  flow 
through  the  tube  must  be  initiated  as  follows: 


a.  Allow  the  apparatus  vapor  temperature  measurement  (channel  40)  to 
reach  at  least  3800  microvolts. 

b.  Cut  in  the  auxiliary  coolant  flow  (50%  or  60%  level)  to  cool  the  vapor 
temperature  to  roughly  3200  microvolts. 

c.  Secure  coolant  flow  through  the  auxiliary  condenser,  and  allow  the 
vapor  temperature  level  to  climb  to  about  3700-3800  microvolts,  which 
allows  a  steam  blanket  to  cover  the  tube. 

d.  Initiate  coolant  flow  through  the  single  horizontal  tube  at  the  80%  level. 
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e.  Cut  in  coolant  flow  to  the  auxiliary  condenser  to  control  pressure,  and 
observe  the  condensation  process  to  ensure  that  a  condensate  film  has 
formed  on  the  tube. 

13.      Run  software  program  DRPINST  for  an  instrumented  tube  (DRPKS  for 
uninstrumented  tubes)  by  pressing  "run"  on  the  keyboard. 

To  take  data  for  an  instrumented  tube,  the  questions  for  DRPINST  can  be 

answered  as  follows: 

•  Select  fluid  ...  Enter  0  for  steam 

•  Select  option  ...  Enter  1  to  take  new  data 

•  Enter  month,  date  and  time  ...  Press  enter 

•  Enter  input  mode  ...  Enter  1  for  new  data 

•  Give  a  name  for  the  raw  data  file  ...  Enter  name 

•  Enter  geometry  code  ...  Enter  1  for  finned,  0  for  plain 

•  Select  insert  type  ...  Enter  0  for  none,  1  for  twisted  tape,  2  for  wire  wrap,  3  for 

Heatex 

•  No.  of  thermocouples  in  wall?  ...  Enter  4,  5,  or  6  depending  on  the  tube 

•  Select  tube  diameter  type  ...  Enter  2  for  medium 

•  Enter  pressure  condition  ...  Enter  0  for  vacuum,  1  for  atmospheric 

•  Give  a  name  for  the  wall  temperature  file  ...  Enter  name 

•  Select  input  ...  Enter  1  for  short,  2  for  long,  or  3  for  raw  data 

•  Like  to  check  NG  (non-condensible  gas)  concentration  ...  Enter  1  for  yes,  2  for 

no;  you  must  answer  yes  for  the  first  data  point 

•  Enter  flowmeter  reading  (%)  ...  Enter  2  digit  number  (i.e.  20  or  58  etc.) 

•  Connect  voltage  line  ...  Flip  the  voltage  line  toggle  switch,  located  on  the  power 

control  panel,  to  the  on  position  and  press  enter 
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•  Disconnect  voltage  line  ...  Flip  the  voltage  line  toggle  switch  off  and  press  enter 

•  Enter  pressure  gauge  reading  (Pga)  ...  Enter  reading  off  gauge  in  psi 

•  Select  measurement  ...  Enter  0  for  teflon,  1  for  metal  sheath,  2  for  quartz, 

3  for  thermopile 

•  Change  TCOOL  rise?  ...  Enter  1  for  yes,  2  for  no 

•  OK  to  store  this  data  set?  ...  Enter  1  for  yes,  0  for  no 

•  Will  there  be  another  run?  ...  Enter  1  for  yes,  0  for  no;  starts  at  check  NG 

concentration  for  following  runs. 

To  take  data  for  an  uninstrumented  tube  the  questions  for  DRPKS  can  be 
answered  as  follows: 

•  Select  fluid  ...  Enter  0  for  steam 

•  Select  option  ...  Enter  0  to  take  new  data 

•  Enter  month,  date  and  time  ...  Press  enter 

•  Enter  disk  number  ...  Enter  number 

•  Enter  input  mode  ...  Enter  0  for  new  data 

•  Select  Ci  ...  Enter  0  to  find  a  Ci  value,  2  to  use  a  Ci  value  stored  in  the 

program. 

•  Give  a  name  for  the  raw  data  file  ...  Enter  name 

•  Enter  geometry  code  ...  Enter  1  for  finned,  0  for  plain 

•  Enter  insert  type  ...  Enter  0  for  none,  1  for  twisted  type,  2  for  wire  wrap,  3  for 

Heatex 

•  Select  tube  type  ...  Enter  0  for  thick  wall  (only  thick  wall  tubes  were  tested) 

•  Select  material  code  ...  Enter  0  for  copper  (only  copper  tubes  were  tested) 
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•  Select  tube  diameter  type  ...  Enter  1  for  medium  (no  small  or  large  diameter 

tubes  were  tested) 

•  Enter  pressure  condition  ...  Enter  0  for  vacuum,  1  for  atmospheric 

•  Want  to  create  a  file  for  NR  vs  F?  ...  Enter  1  for  yes,  0  for  no 

•  Give  a  name  for  plot  data  file  ...  Enter  name;  easiest  to  use  the  raw  data  file 

name  preceded  by  a  P 

•  Select  output ...  Enter  0  for  short,  1  for  long,  2  for  raw  data 

•  Like  to  check  NG  concentration  ...  Enter  1  for  yes,  2  for  no;  you  must  answer 

yes  for  the  first  data  point 

•  Enter  flowmeter  reading  (%)  ...  Enter  2  digit  number  (i.e.  20  or  60  etc.) 

•  Connect  voltage  line  ...  Flip  the  voltage  line  toggle  switch  on  and  press  enter 

•  Disconnect  voltage  line  ...  Flip  the  voltage  line  toggle  switch  off  and  press  enter 

•  Enter  pressure  gauge  reading  (Pga)  ...  Enter  reading  off  gauge  in  psi 

•  Select  measurement  ...  Enter  0  for  teflon,  1  for  metal  sheath,  2  for  quartz, 

3  for  thermopile 

•  Change  TCOOL  rise  ...  Enter  1  for  yes,  2  for  no 

•  OK  to  store  this  data  set  ...  Enter  1  for  yes,  0  for  no 

•  Will  there  be  another  run  ...  Enter  1  for  yes,  0  for  no;  starts  at  check  NG 

concentration  for  following  runs. 


14.  Only  answer  the  program  questions  up  to  "Enter  flowmeter  readings". 
Monitor  system  temperature  using  the  vapor  thermocouple  voltage  reading 
(the  program  automatically  resets  to  channel  40)  closely  until  system  warmup 
is  complete. 

15.  Monitor  system  temperature  and  pressure  carefully  to  prevent  a  system  over 
pressure  during  warmup  (especially  at  atmospheric  conditions). 
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16.  If  conducting  a  vacuum  run,  gradually  adjust  voltage  to  90  volts  (usually  in  10 
volt  increments).  Obtain  the  desired  operating  condition  by  manually 
controlling  coolant  flow  through  the  auxiliary  condenser  until  channel  40 
reads  1970  ±  20  microvolts  (— 48°C).   Vapor  velocity  ~  2  m/s. 

17.  If  conducting  an  atmospheric  run,  gradually  adjust  voltage  to  175  volts  from 
the  90  volt  level  in  10-20  volt  increments.  Again  the  desired  operating 
condition  is  obtained  by  manually  controlling  coolant  flow  through  the 
auxiliary  condenser  until  channel  40  reads  4280  ±  20  microvolts  (~100°C). 
Vapor  velocity  ~  1  m/s. 

18.  Monitor  the  condensation  process  using  the  glass  viewing  window  periodically 
to  ensure  that  filmwise  condensation  is  maintained.  To  clear  the  viewing 
window  of  fog  and  moisture  increase  coolant  flow  through  the  auxiliary 
condenser  briefly  to  50%  or  60%,  then  reset  to  desired  flow  rate. 

19.  When  taking  readings  be  sure  to  check  the  flowmeter  setting  prior  to  entering 
it  into  the  computer  (it  has  a  tendency  to  fluctuate  slightly). 

20.  If  conducting  vacuum  and  atmospheric  runs  on  the  same  day  always  conduct 
the  vacuum  run  first.  If  the  atmospheric  run  is  done  first  it  takes  too  long  for 
the  system  to  cool  down  to  vacuum  operating  temperatures. 


The  system  is  secured  in  the  following  manner: 


1.  Secure  power  to  the  heating  elements. 

2.  Secure  coolant  flow  through  the  tube,  through  the  auxiliary  condenser,  and 
to  the  sump  tank. 

3.  If  desired  to  maintain  the  system  at  vacuum  conditions  until  the  next  run  the 
shutdown  is  complete.  Continued  cooling  water  circulation  may  be  used  to 
assist  in  cooling  down  the  system. 

4.  To  bring  the  system  back  to  atmospheric  conditions  slowly  open  the  vent 
valve. 

5.  The  data  acquisition  system  may  be  turned  off  whenever  it  is  not  necessary 
to  monitor  system  parameters. 

6.  Periodically  change  distilled  water  in  the  boiler. 
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7.  If  an  emergency  should  arise  such  as  abrupt  overpressurization  or  breakage, 
immediately  secure  power  to  the  heaters  and  open  the  vent  valve,  then  let  the 
system  cool  down  before  checking  the  apparatus  for  damage. 

B.       EXPERIMENTAL  PROCEDURES  AND  OBSERVATIONS 

Water  is  a  poor  wetting  medium  and  therefore  great  care  was  taken  to  ensure 
that  uniform  filmwise  condensation  was  the  only  condensation  mode  occurring  during 
a  data  run.  Even  though  the  apparatus  was  meticulously  cleaned  (as  mentioned 
previously),  a  continuing  problem  with  dropwise  condensation  manifested  itself. 
Subsequent  to  steam  cleaning  the  system  with  a  Sparkleen  soap  solution,  by 
operating  the  system  with  a  soapy  solution  in  the  boiler  (the  solution  bubbled 
through  the  entire  apparatus),  dropwise  condensation  was  observed  on  the  installed 
instrumented  tube.  After  taking  some  data  when  in  the  dropwise  condition,  the  tube 
was  removed  and  rigorously  cleaned  using  a  warm  Alconox  soap  solution  with  a 
scrub  brush.  However,  after  observing  the  filmwise  mode  initially,  the  condensation 
mechanism  soon  transitioned  to  mixed  mode  and  then  back  to  the  dropwise  mode. 

Since  only  a  filmwise  condition  over  several  hours  would  suffice,  the  tube 
chemical  treatment  procedure  used  by  Guttendorf  [Ref.  3]  and  several  other 
researchers  was  used  to  produce  filmwise  condensation.  The  tube  was  chemically 
treated  prior  to  installation  as  follows: 


1.  Clean  the  internal  and  especially  the  external  surfaces  of  the  tube  using  a  soft 
brush  and  mild  soap  (using  the  Alconox  detergent  in  warm  water),  rinse  with 
acetone  then  rinse  thoroughly  with  distilled  water.    Repeat  the  cleaning 
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procedure  until  the  distilled  water  rinse  perfectly  wets  the  tube  surface;  any 
breaks  in  the  wetting  film  at  this  point  are  likely  to  result  in  dropwise 
condensation  spots  once  the  tube  is  installed  in  the  apparatus. 

2.  Place  the  tube  in  a  steam  bath. 

3.  Mix  equal  amounts  of  ethyl  alcohol  and  a  50%  by  weight  solution  of  sodium 
hydroxide.  Keep  the  solution  warm  so  that  a  watery  consistency  is 
maintained. 

4.  Apply  the  solution  to  the  tube  with  a  small  paint  brush,  retaining  the  tube  in 
the  steam  bath.  If  the  tube  has  not  been  treated  previously,  apply  a  coating 
of  the  solution  every  10  minutes  for  an  hour.  If  the  tube  has  been  previously 
treated,  apply  a  coating  every  5  minutes  for  a  period  of  20  minutes. 

5.  Remove  the  tube  from  the  steam  bath  and  thoroughly  rinse  the  tube  with 
distilled  water  to  remove  any  excess  solution.  Install  the  tube  in  the  test 
section  immediately,  being  careful  not  to  touch  the  tube  surface.  Oil  or  dirt 
from  any  source  may  contaminate  the  tube  surface  and  result  in  mixed  mode 
or  dropwise  condensation. 


The  oxide  layer  which  forms  on  the  tube  is  very  thin,  and  has  negligible  thermal 
resistance  and  high  wetting  characteristics. 

Once  the  tube  was  installed  in  the  apparatus  (with  the  desired  insert  in  place), 
the  system  startup  procedure  outlined  in  section  IV  A  was  followed  to  take  data  at 
desired  conditions. 

At  vacuum  conditions,  when  single  tube  coolant  flow  was  initiated  with  vapor 
velocity  at  ~  2  m/s,  the  condensation  on  the  tube  did  not  develop  as  a  perfect  film 
but  instead  left  patches  where  the  film  was  broken.  These  patches,  or  streaks 
seemed  to  occur  at  regular  intervals,  and  it  was  postulated  that  they  were  due  to 
vortex  shedding  of  vapor  around  the  tube.  Therefore,  the  procedure  in  the  startup 
section  IV  A,  step  12,  was  used  to  promote  the  development  of  uniform  filmwise 
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condensation  by  inducing  a  stationary  vapor  condition  around  the  tube.  This  allowed 
a  steam  blanket  to  form  around  the  tube  prior  to  coolant  flow  initiation.  After  flow 
initiation  the  appearance  of  the  condensate  film  on  the  tube  surface  was  continuous 
with  no  breaks.  Momentary  film  instabilities  were  observed  at  vacuum  conditions  at 
pressures  below  ■—  20  kPa  after  a  continuous  film  was  established  at  higher  pressures. 
These  instabilities  seemed  to  interrupt  the  film  sheet  only  for  an  instant  and  then 
disappear,  and  may  have  been  caused  by  vortex  shedding  of  the  vapor  around  the 
tube  as  already  mentioned.  A  possible  mechanism  to  explain  these  instabilities  is 
that  the  higher  vapor  velocity  at  vacuum  conditions  momentarily  thins  the  condensate 
film  via  vortex  shedding,  yet  this  thinning  effect  is  overcome  by  surface  tension  forces 
in  the  film  sheet  which  tend  to  restore  the  continuous  film.  These  instabilities  could 
only  be  seen  for  an  instant  and  then  would  vanish,  being  very  transitory  in  nature. 
There  also  did  not  appear  to  be  any  pattern  whatsoever  to  the  instability  formation. 
The  instabilities  were  not  observed  at  pressures  above  —  20  kPa. 

The  data  taking  regimen  for  each  data  set  involved  starting,  then  verifying  the 
existence  of  a  filmwise  condensation  condition,  then  taking  data  at  flow  rates  (in  %) 
of  80,  70,  60,  50,  40,  30,  and  20  then  back  to  80  and  50  to  check  for  repeatability 
within  the  data  set.  Two  data  points  were  taken  at  each  of  the  first  seven  data  points 
and  one  each  for  the  last  two,  which  gave  a  total  of  16  data  points.  It  was  usually 
quite  clear  from  the  two  comparison  points,  and  from  data  taken  previously  under 
similar  conditions,  whether  the  data  set  should  be  rejected  or  accepted.  After  tube 
installation,  the  appearance  of  one  or  more  small  patches  (breaks  in  the  film)  after 
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—  7-10  hours  of  operation  signaled  the  beginning  of  tube  contamination  which  got 
worse  with  time.   The  tube  would  then  be  removed  and  cleaned. 

C.      TUBES  TESTED 

The  data  taken  during  this  study  involved  extensive  use  of  the  instrumented 
smooth  tube  (SOI)  fabricated  by  Poole  [Ref.  6],  (six  wall  thermocouples  spaced  60 
degrees  apart  placed  at  midwall  and  midlength).  Due  to  excessive  thermocouple 
wear,  only  5  thermocouples  in  this  instrumented  tube  functioned  properly.  The  tube 
was  positioned  in  the  apparatus  to  make  optimum  use  of  functioning  thermocouples. 
The  preferred  arrangement  placed  4  thermocouples  at  10°,  190°,  250°,  and  310°  from 
the  top  dead  center  position  of  the  tube.  This  arrangement  provided  readings  from 
the  top  and  bottom  of  the  tube  and  two  intermediate  points,  giving  the  most  accurate 
mean  tube  wall  temperature  and  best  temperature  profile  readings.  At  the 
conclusion  of  this  study,  the  manufacture  of  the  new  instrumented  tubes  was  not  yet 
complete;  progress  to  date  (March  1992)  is  recorded  in  Appendix  F. 

Data  was  also  taken  on  a  uninstrumented  smooth  tube  (S02)  and  four  finned 
tubes  (S03,  S04,  S05,  S06)  with  fin  spacings  of  0.5  mm,  1.0  mm,  and  1.5  mm,  and  2.0 
mm.  All  tubes  tested  were  classified  as  medium  tubes,  with  an  outside  diameter  of 
19.05  mm,  and  an  inside  diameter  of  12.70  mm.  The  four  finned  tubes  all  had  the 
same  fin  height  of  1.0  mm,  and  fin  width  of  1.0  mm.  Data  runs  were  taken  either 
with  no  insert,  with  the  wire  wrap  or  the  Heatex  insert  installed.By  more  efficiently 
mixing  the  coolant  an  insert  significantly  increases  the  inside  heat-transfer  coefficient. 
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The  Heatex  insert  consists  of  a  central  wire  core  onto  which  are  wound  a  series 
of  wire  loops,  each  inclined  at  a  common  angle  to  the  core.  The  loops  come  into 
direct  contact  with  the  tube  wall,  and  each  loop  provides  a  significant  amount  of 
coolant  mixing  as  the  coolant  flows  through  the  loop  mesh.  [Ref.  12]. 

The  wire  wrap  insert  was  a  copper  wire  spirally  wrapped  around  a  central 
stainless  steel  rod  with  a  uniform  pitch.  This  insert  induced  a  swirling  coolant 
motion,  which  enhanced  turbulent  mixing  within  the  tube.  This  particular  wire  wrap 
insert  was  the  same  insert  used  by  Guttendorf  [Ref.  3],  Coumbes  [Ref.  22],  and  Van 
Petten  [Ref.  2]  to  collect  data  on  the  medium  family  of  tubes. 

A  summary  of  data  runs  is  given  in  Table  2. 

Table  2.   DATA  RUN  SUMMARY 


Tube 

Filename 

Pressure 
(kPa) 

Vapor 
Velocity 

(m/s) 

Insert 

SOI 

DSOIMVSH2 

11 

2.17 

Heatex 

DSOIMASH4 

101 

1.06 

Heatex 

DSOIMASH5 

101 

1.06 

Heatex 

FSOIMASH3 

11 

2.11 

Heatex 

FIMAVSH1 

28 

3.53 

Heatex 

FIMAVSH2 

41 

2.45 

Heatex 

FIMAVSH3 

69 

1.53 

Heatex 

FIMAVSH4 

101 

1.08 

Heatex 

FIMAVSH5 

28 

1.01 

Heatex 

FIMAVSH6 

41 

1.04 

Heatex 

FIMAVSH7 

69 

1.02 

Heatex 

FIMASW3 

101 

1.07 

wire  wrap 

FIMASW4 

101 

1.07 

wire  wrap 
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Tube 

Filename 

Pressure 
(kPa) 

Vapor 
Velocity 

(m/s) 

Insert 

FIMASW5 

101 

1.06 

wire  wrap 

COMPAIW1 

101 

1.07 

wire  wrap 

COMPAIW2 

101 

1.07 

wire  wrap 

FIMASN4 

101 

1.08 

none 

FIMASN5 

101 

1.07 

none 

FIMASN6 

101 

1.08 

none 

COMPAIN1 

101 

1.08 

none 

COMPAIN2 

101 

1.08 

none 

S02 

FNMAVSH1 

15 

6.24 

Heatex 

FNMAVSH2 

41 

2.42 

Heatex 

FNMAVSH3 

69 

1.51 

Heatex 

FNMAVSH4 

102 

1.03 

Heatex 

FNMAVSH5 

28 

1.02 

Heatex 

FNMAVSH6 

41 

1.03 

Heatex 

FNMAVSH7 

69 

1.02 

Heatex 

FNMAVSH8 

101 

1.07 

Heatex 

S03 

FIMAF051 

101 

1.06 

Heatex 

FIMAF052 

101 

1.06 

Heatex 

S04 

FIMAF101 

101 

1.07 

Heatex 

FIMAF102 

101 

1.06 

Heatex 

S05 

FIMAF151 

101 

1.07 

Heatex 

FIMAF152 

101 

1.07 

Heatex 

S06 

FIMAF201 

101 

1.07 

Heatex 

FIMAF202 

101 

1.07 

Heatex 
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Tube  descriptions: 

501  instrumented  smooth  tube  fabricated  by  Poole 

502  uninstrumented  smooth  tube 

503  0.5  mm  finned  tube 

504  1.0  mm  finned  tube 

505  1.5  mm  finned  tube  (this  tube  designated  as  F096  by  Guttendorf;  also  a 
similar  1.5  mm  finned  tube  designated  as  F006  by  Guttendorf  and  Van 
Petten  was  not  tested);  the  1.5  mm  finned  tube  results  were  not  included 
in  subsequent  finned  tube  analysis  since  the  reason  for  differences  in 
experimental  results  for  these  two  tubes  has  not  yet  been  resolved. 

506  2.0  mm  finned  tube 
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V.   THEORETICAL  BACKGROUND  AND  DATA  REDUCTION  PROCEDURES 

A.       THEORETICAL  BACKGROUND 

The  overall  or  total  resistance  to  heat  transfer  from  vapor  to  coolant  consists 
of  the  sum  of  the  vapor  side  resistance  (R0),  the  tube  wall  resistance  (Rv,),  and  the 
coolant  side  resistance  (R,);  this  neglects  any  fouling  resistance  since  clean  tubes  are 
always  used. 

R     ,  =  R  +  R  +R  (5.1) 

total        o        w        i  v       ' 

The  vapor  and  coolant  side  resistances  are  convective  in  nature  and  may  be 
expressed  by  the  reciprocal  of  their  respective  heat-transfer  coefficient  and  surface 
area  product. 


(5.2) 


K- 

o     0 

K- 

1 

(5.3) 


i    < 


where: 

R0  =  outside  vapor  side  resistance  to  heat  transfer  (K/W) 

h0  =  outside  heat-transfer  coefficient  (W/m2K) 

A0  =  effective  outside  surface  area  (m  ) 

R,  =  inside  coolant  side  resistance  to  heat  transfer  (K/W) 

h,  =  inside  heat-transfer  coefficient  (W/m2K) 

A,  =  effective  inside  surface  area  (m  ) 

The  tube  wall  resistance  is  conductive  in  nature  and  is  represented  by  the  radial 

conduction  equation. 


36 


■3 

(5.4) 


W 


2nLk 


where: 


R^,  =  tube  wall  resistance  (K/W) 

Dr  =  outside  or  root  (if  finned)  diameter  (m) 

D,  =  inside  tube  diameter  (m) 

L  =  active  condensing  length  (133  mm) 

1^  =  thermal  conductivity  of  tube  wall  (W/m-K) 

The  effective  outside  area  of  the  tube  is  calculated  using  the  following 

expression: 

Ao  =  nDrL  (5.5) 

where: 

Dr      =    outside  or  root  (if  finned)  diameter  (m) 
L        =    active  condensing  length  (133  mm) 

The  effective  inside  area  includes  the  inside  surface  area  involving  the  active 

condensing  length  and  the  inside  surface  area  of  the  insulated  inlet  and  outlet 

portions  of  the  tube.    These  portions  of  the  tube  act  as  fins  and  remove  heat  via 

axial  conduction.    The  extended  fin  assumption  with  associated  fin  efficiencies  was 

used  to  account  for  these  end  losses. 


^  =  7iZ).(L  +  lllL1+Ti2L2)  (5.6) 


where: 


D,  =  inside  diameter  of  tube  (m) 

L  =  active  tube  condensing  length  (m) 

17!  =  fin  efficiency  of  inlet  portion  of  tube 

Lx  =  length  of  inlet  portion  of  tube  (m) 

172  =  fin  efficiency  or  outlet  portion  of  tube 

L2  =  length  of  outlet  portion  of  tube  (m) 

The  overall  thermal  resistance  to  heat  transfer  may  be  expressed  by: 
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1 


UA, 


(5.7) 


o     o 


Substituting  equations  (5.2)  and  (5.3)  into  equation  (5.7)  yields: 


1  1         D  1 

UA      hA^      w    hA 


(5.8) 


o     o 


o     o 


I     I 


where: 

U0  =    overall  heat-transfer  coefficient  (W/m2K) 

A0  =    effective  outside  surface  area  (m  ) 

h0  =    outside  heat-transfer  coefficient  (W/m2K) 

R^  =    Tube  wall  resistance  (K/W) 

h,  =    inside  heat-transfer  coefficient  (W/m2K) 

The  single  tube  condenser  apparatus  uses  the  log  mean  temperature  difference 

(LMTD)  analysis  for  calculation  of  the  heat  transfer  between  the  hot  vapor  and  cold 

coolant. 


Q-UAALMTD) 


o     o 


(5.9) 


where: 


Q 

A„ 


=  heat  transfer  rate  to  the  cooling  water 
=  overall  heat-transfer  coefficient  (W/m2 
=    effective  outside  surface  area  (m  ) 


•S? 


LMTD    =    log  mean  temperature  difference  between  vapor  and  coolant  (K) 
The  log  mean  temperature  difference  (LMTD)  is  given  by: 


LMTD  = 


In 


T    -T 


(5.10) 
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where: 

Tt       =    coolant  inlet  temperature  (K) 
T2      =    coolant  outlet  temperature  (K) 
Tsat     =    vapor  saturation  temperature  (K) 

In  this  and  previous  studies  at  NPS,  the  quartz  thermometer  output  for  Tt  and 
T2  were  used  in  the  calculations  for  the  coolant  temperature  rise,  and  the  saturation 
temperature,  Tsat,  was  measured  using  the  vapor  thermocouple  (channel  40). 

The  total  heat  transfer  across  the  tube  is  experimentally  determined  by 
measuring  the  mass  flow  rate  of  fluid  through  the  tube  and  its  accompanying 
temperature  rise. 

Q  =  mc(T2-Tl)  (5.11) 


where: 

Q  =  heat  transfer  rate  (W) 

m  =  mass  flow  rate  of  coolant  (kg/s) 

Cp  =  specific  heat  of  coolant  at  constant  pressure  (J/kg-K) 

Tx  =  coolant  inlet  temperature  (K) 

T2  =  coolant  outlet  temperature  (K) 

Equation  (5.11)  may  be  used  directly  from  the  experimental  data.     The 

resultant  heat  transfer  rate,  Q,  is  then  substituted  into  equation  (5.9)  to  find  the 

overall  heat-transfer  coefficient,  U0. 

U=   ,    ,.?__.  (5.12) 


A0(LMTD) 


where: 


Q       =    heat  transfer  rate  from  eq.  (5.11)  (W) 

A0      =    effective  outside  surface  area  (m2) 

LMTD    =    log  mean  temperature  difference;  eq.  (5.10)  (K) 
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Since  R^  U0,  A,,,  and  A,  are  known  quantities,  this  leaves  only  two  unknowns 
in  equation  (5.8),  h0  and  hy,  the  outside  and  inside  heat  transfer  coefficients. 

Often  the  coolant  side  thermal  resistance  is  dominant,  and  inserts  such  as  those 
mentioned  in  section  IV  C  are  used  to  lower  the  inside  resistance.  This  allows  a 
more  accurate  computation  of  the  outside  heat-transfer  coefficient,  h0,  when  using 
the  modified  Wilson  plot  technique  mentioned  in  section  V  B.  Vapor  side  heat 
transfer  may  also  be  enhanced  through  the  use  of  fins,  drainage  strips,  or  dropwise 
condensation  promoters. 

B.       MODIFIED  WILSON  PLOT  TECHNIQUE 

The  ideal  way  to  solve  for  h0  and  h,  in  equation  (5.8)  is  through  the  use  of 
instrumented,  which  accurately  determine  a  mean  tube  wall  temperature.  The  inside 
and  outside  mean  tube  wall  temperatures  may  then  be  obtained  directly  by  assuming 
a  linear  temperature  profile  across  the  wall.  Since  the  vapor  temperature  and  mean 
coolant  temperature  are  known,  the  inside  and  outside  heat-transfer  coefficients  may 
then  be  calculated  directly  using  equation  (5.13). 

q  =  hAT  (513) 

where: 

q        =    heat  flux  (W/m2) 

h        =    heat-transfer  coefficient  (h,  for  inside,  h0  for  outside)  (W/m2K) 
AT     =    temperature  difference  across  resistive  medium  (AT  =  Tsat-Twanoutside 
for  the  vapor  side,  and  AT  =  TwaUinside  -  T^,,,  for  the  coolant  side) 

For  data  collection  on  a  large  number  of  tubes,  the  use  of  instrumented  tubes 
is  impractical  due  to  the  high  cost  and  difficulty  involved  in  manufacturing  so  many 
tubes.  Therefore,  the  modified  Wilson  plot  technique  was  developed,  which  solves 
for  the  inside  and  outside  heat-transfer  coefficients  simultaneously  without  using  wall 
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thermocouples.  To  obtain  the  most  accurate  results  with  this  method,  it  is  necessary 
that  the  inside  and  outside  coefficients  be  relatively  equal  in  magnitude. 

The  modified  Wilson  plot  technique  requires  that  the  "form"  of  the  equation 
for  both  the  inside  and  outside  heat-transfer  coefficients  be  known.  The  Nusselt 
theory  and  Sieder-Tate  correlation  are  used  to  represent  the  "form"  of  the  outside 
and  inside  heat-transfer  coefficients  respectively.  The  Nusselt  theory  when  based  on 
q  can  be  represented  by: 


ho  =  a 


kjpjgh 


fg 


VfVoV 


1/3 


=aZ 


(5.14) 


or 


where: 


a 
Pf 

q 


=  outside  heat-transfer  coefficient  (W/m2K);  based  on  q 

=  dimensionless  coefficient 

=  thermal  conductivity  of  condensate  film  (W/m-K) 

=  condensate  film  density  (kg/m3) 

=  gravitational  constant  (9.81  m/s2) 

=  specific  enthalpy  of  vaporization  (J/kg) 

=  dynamic  viscosity  of  condensate  film  (Ns/m2) 

=  outside  tube  diameter  (m) 

=  heat  flux  based  on  outside  area  (Q/A0)  (W/m2) 


The  Sieder-Tate  correlation  may  be  represented  by: 


h=C-^Re°*Pr113 


/\LC\0U 


i         i 


Z) 


IP 


=  Cfi 


w/ 


(5.15) 


where: 


Q 

kc 

D, 

Re 

Pr 


=  inside  heat  transfer  coefficient  (W/m2K) 

=  Sieder-Tate  leading  coefficient 

=  thermal  conductivity  of  cooling  water  (W/m-K) 

=  inside  tube  diameter  (m) 

=  Reynolds  number 

=  Prandtl  number 

=  dynamic  viscosity  of  cooling  water  at  bulk  temperature  (Ns/m2) 
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fiw      =    dynamic  viscosity  of  cooling  water  at  mean  inner  wall  temperature 
(N-s/m2) 

Substituting  equations  (5.14)  and  (5.15)  into  equation  (5.8)  gives: 


U. 


M 


AnZ 


1 

+  — 
CiQAi     a 


(5.16) 


By  letting: 


Y  = 


1 

-K*. 

x= 

(5.17) 


(5.18) 


m  = 


(5.19) 


b-1 

a 


(5.20) 


then  equation  (5.16)  reduces  to: 


Y=mX+b 


(5.21) 


The  parameters  CI  and  Z  are  both  temperature  dependent,  and  must  be 
determined  iteratively.  A  least-squares  fit  of  equation  (5.21)  is  used  to  determine 
C,  and  a.  Once  C,  is  known,  h,  can  be  calculated  using  equation  (5.15).  With  h,  and 


42 


U0  known,  the  value  of  h0  can  then  be  easily  determined  by  rearranging  equation 
(5.7)  as  in  equation  (5.22),  or  by  using  equation  (5.14)  with  the  value  of  a  known. 


J_     J_ 

o  o 


Mi 


°-+K\ 


(5.22) 


C.       INSTRUMENTED  TUBE  IMPROVEMENTS  FOR  DATA  REDUCTION 

In  previous  work  at  NPS,  the  standard  form  of  the  Sieder-Tate  equation  was 
used  with  a  Reynolds  number  exponent  of  0.8,  as  in  equation  (2.10).  One  of  the 
aims  of  this  thesis  was  to  use  an  instrumented  tube  to  directly  determine  the  inside 
and  outside  coefficients,  h,  and  h0  ,  and  then  use  the  data  to  determine  a  more 
"exact"  form  of  the  Sieder-Tate-type  equation  to  be  used  for  each  insert.  The 
coolant  side  correlations  mentioned  in  section  II  C  were  based  on  a  long,  straight 
entrance  length.  The  sharp  90°  bend  just  prior  to  the  test  section  tube  entrance 
undoubtedly  creates  entrance  effects  which  lead  to  discrepancies  between  our 
experimental  data  and  heat  transfer  behavior  predicted  by  the  Sieder-Tate 
correlation. 

Assuming  the  final  form  of  the  inside  heat-transfer  correlation  to  be  that  of 
equation  (2.10)  gives  the  following: 


Nu=CRemPrlf3 


(Vet1* 


^w 


(5.23) 


where: 

m       =    Reynolds  number  exponent  to  be  determined 
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Rearranging  equation  (5.23)  gives: 


Nu 


PrW\ 


(\lc\014 


=  C:Rem 


wj 


(5.24) 


Taking  the  natural  log  of  equation  (5.24)  yields: 


In 


Nu 


jyi/3 


-  In  Ci+  m  In  Re 


(5.25) 


Equation  (5.25)  is  in  the  form  of  a  linear  equation,  and  by  plotting  In 
(Nu/(PTm(jiJ/jiw)ou))  versus  In  Re,  the  slope  and  intercept,  namely  the  Reynolds 
number  exponent  and  Sieder-Tate  coefficient  may  be  determined  from  the 
instrumented  data. 

With  the  unknown  parameters  of  equation  (5.23)  determined,  the  new  inside 
heat-transfer  correlations  (one  for  each  insert)  could  be  used  in  the  data  reduction 
program  to  give  the  value  of  the  inside  heat- transfer  coefficient,  h„  directly  and  to 
provide  a  more  accurate  calculation  of  the  outside  heat-transfer  coefficient  h0. 

D.      ENHANCEMENT  RATIO 

Following  the  development  of  Van  Petten  [Ref.2],  and  Nusselt  theory, 
experimental  data  can  be  curve  fitted,  using  a  least-squares  analysis,  to  an  equation 
of  the  following  form: 
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q  =  a&Tn 


(5.26) 


where: 


q        =    heat  flux  based  on  outside  area  (Q/A0)    fW/m2) 
AT     =    temperature  drop  across  the  condensate  film  (K) 

Substituting  equation  (5.12)  into  his  expression  yields: 


h    =  aLT 

o 


n-\ 


(5.27) 


From  Nusselt  theory  n  =  0.75,  therefore  the  enhancement  ratio,  based  on 
constant  temperature  drop  across  the  condensate  film,  can  be  expressed  as: 


h4   -    af 


AT 


(5.28) 


OS 


a. 


where: 

6 


AT 


/ 
S 

a 


=    enhancement  ratio  based  on  constant  temperature  drop  across  the 

condensate  film 
=    subscript  denoting  finned  tube 
=    subscript  denoting  smooth  tube 
=    outside  heat-transfer  coefficient  (W/m2K) 
=    constant  of  proportionality  introduced  in  equation  (5.26) 


Also  for  constant  AT;  using  equation  (5.14): 


s  s 


*/ 


1/3    aA 


a. 


*A 


J]  1/3 


a/iT 


.  af 


a. 


a. 


7J 


1/3 


(5.29) 


If  the  heat  flux  is  kept  constant  the  values  of  Z/  and  Zs  remain  equal,  which 
results  in  equation  (5.30). 
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e  =  !k\  -  afzf  -  af 

q  hos   q      asZs       as 


(5.30) 


Combining  equations  (5.29  and  (5.30)  gives  the  relationship  between  eAT  and 
€   in  equation  (5.33). 


ff 

S 


a. 


a, 


W_af 
a. 


(5.31) 


a 

s 


a. 


w 


(5.32) 


€    =(€      )4/3 


(5.33) 


Note  that  eAT  and  e   are  independent  of  q  and  AT. 
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VI.    RESULTS  AND  DISCUSSION 

A.       DROPWISE  CONDENSATION 

As  mentioned  in  section  IV  B,  a  dropwise  condensation  condition  was  obtained 
initially.  Data  were  taken  during  dropwise  conditions  at  vacuum  and  atmospheric 
pressure  with  the  instrumented  tube,  and  this  data  is  compared  to  filmwise  data  in 
Figure  4.  The  figure  shows  a  marked  contrast  between  filmwise  and  dropwise 
condensation  data.  The  dropwise  heat-transfer  enhancement,  compared  to  the 
filmwise  data,  varied  from  —2  to  7  for  vacuum  conditions  and  —9  to  10  for 
atmospheric  conditions.  Marto  et  al  [Ref.  23]  studied  the  use  of  organic  coatings  for 
the  promotion  of  dropwise  condensation  of  steam,  and  obtained  an  outside  heat- 
transfer  coefficient  of  —55  k\V/m2K  for  a  Fluoroacrylic  coating  compared  to  —30  to 
85  kW/m2K  found  under  vacuum  conditions  in  this  study  (both  taken  at  P—  11  kPa 
and  vapor  velocity  —2  m/s).  The  difficulty  in  accurately  measuring  the  temperature 
drop  across  the  condensate  film  for  the  dropwise  condition  is  illustrated  by  the  large 
amount  of  scatter  in  the  dropwise  data;  this  problem  is  caused  by  the  lack  of  a  stable 
film  and  the  intermittent  presence  of  drops  near  the  instrumented  tube  wall 
thermocouples.  The  exact  cause  of  the  dropwise  conditions  was  never  determined. 
Clearly,  some  organic  contamination  either  from  the  boiler  feed  water  or  from  the 
gasket  material  was  depositing  on  the  test  tube.  Since  this  thesis  was  devoted  to 
filmwise  condensation,  great  efforts  were  made  to  clean  the  test  tube  and  prepare 
it  chemically  so  that  the  condensate  would  wet  the  surface. 
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Figure  4.  Comparison  of  Dropwise  and  Filmwise  Condensation  Data  (Smooth 
Instrumented  Tube,  Heatex  Insert) 
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B.       INSTRUMENTED  TUBE  RESULTS 

After  resorting  to  the  oxide  coating  procedure  mentioned  previously  to  mitigate 
the  dropwise  contamination  problem,  a  series  of  runs  were  made  with  the 
instrumented  tube  fabricated  by  Poole  [Ref.  6]  with  the  Heatex  insert  installed. 
Using  a  mean  wall  temperature,  the  inside  and  outside  heat-transfer  coefficients 
could  be  evaluated  directly. 

Figure  5  shows  filmwise  condensation  data  taken  at  various  pressures  and  vapor 
velocities.  The  four  data  sets  for  increasing  vapor  velocity  (from  1  to  3.5  m/s)  and 
decreasing  vapor  pressure  (from  101  to  28  kPa)  show  the  effect  of  vapor  shear 
thinning  the  condensate  film,  giving  an  increase  in  the  outside  heat-transfer 
coefficient.  This  series  of  data  runs  was  taken  by  maintaining  the  heater  voltage  at 
175  volts  and  adjusting  coolant  flow  through  the  auxiliary  condenser  to  control 
pressure. 

The  four  data  sets  in  Figure  5  with  constant  vapor  velocity  (—1  m/s)  and 
decreasing  vapor  pressure  (from  101  to  28  kPa)  show  reduction  of  the  outside  heat- 
transfer  coefficient  with  decreasing  saturation  temperature  (due  to  decreasing 
saturation  pressure).  This  effect  is  thought  to  be  due  to  an  increase  in  condensate 
viscosity  at  lower  temperatures  which  tends  to  prevent  the  condensate  from  flowing 
around  and  draining  from  the  tube  as  easily  as  at  higher  temperatures.  The  resulting 
condensate  film  thickening  provides  an  additional  resistance  to  heat  transfer  thereby 
lowering  the  outside  heat-transfer  coefficient.  This  series  of  data  runs  was  taken  by 
both  adjusting  heater  power  and  auxiliary  condenser  coolant  flow  to  obtain  the  same 
pressure  conditions  as  above  with  a  constant  vapor  velocity  of  —I  m/s. 
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Figure  5.   Effect  of  Pressure  and  Vapor  Velocity  on  the  Steam 

Heat-Transfer  Coefficient  (Smooth  Instrumented  Tube, 
Heatex  Insert) 
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Figure  6  shows  how  increased  pressure,  with  vapor  velocity  held  constant,  gives 
a  larger  increase  in  the  outside  heat-transfer  coefficient  over  the  predicted  Nusselt 
value.  This  discrepancy  is  also  thought  to  be  caused  by  the  viscosity  effect  since  the 
Nusselt  treatment  seemingly  takes  into  account  the  other  possible  causes.  Also 
depicted  more  clearly  is  the  dramatic  enhancement  due  to  vapor  shear  effects  for 
increased  vapor  velocity  at  constant  pressure. 

As  mentioned  previously,  the  use  of  inserts  allows  more  effective  mixing  of 
coolant  and  facilitates  greater  accuracy  in  the  calculation  of  the  outside  heat-transfer 
coefficient.  Figure  7  shows  the  outside  heat-transfer  coefficient  determined  from  the 
use  of  the  two  inserts  together  with  no  insert.  The  results  for  the  wire  wrap  and 
Heatex  inserts  are  closely  grouped,  whereas  the  no  insert  data  shows  much  greater 
scatter.  This  seems  to  indicate  that  the  use  of  an  effective  insert  does  indeed 
enhance  accuracy  for  the  instrumented  tube  data. 

The  mean  temperature  difference  across  the  condensate  film  for  the  no  insert 
case  has  a  significantly  lower  value  than  either  of  the  insert  cases;  with  the  outside 
and  tube  wall  conditions  unchanged,  this  indicates  a  higher  inside  resistance  for  the 
no  insert  case.  To  illustrate  this  point  a  mid-range  instrumented  tube  data  point  at 
the  same  conditions  for  each  insert  case  is  shown  in  Table  3  (P^,— 101  kPa, 
Vvapor~  1.1  m/s,  Vcoolant~2.75  m/s).  It  shows  that  both  inserts  have  a  comparable  effect 
(slightly  better  with  Heatex)  and  roughly  provide  a  factor  of  two  enhancement  over 
the  no  insert  case  in  the  inside  heat  transfer  coefficient.  Due  to  the  increased  inside 
resistance  for  the  no  insert  case,  the  heat  flux  shows  a  decrease  of  about  15%. 
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Figure  6.   Comparison  of  Experimental  Results  with  Nusselt  Theory 
for  Varying  Pressures  and  Vapor  Velocities  (Smooth 
Instrumented  tube,  Heatex  Insert) 
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Figure  7.   Comparison  of  Steam  Heat-Transfer  Coefficients  at  Atmospheric 

Conditions  (P^—lOl  kPa,  Vvapor~l.l  m/s)  for  Three  Insert  Conditions 
(Smooth  Instrumented  Tube) 
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Table  3.  PERFORMANCE  COMPARISON  OF  TUBE  INSERTS 
AT  ATMOSPHERIC  CONDITIONS 


Heatex  Insert 

Wire  Wrap  Insert 

No  Insert 

h,  (kW/m2K) 

30.85 

27.63 

15.31 

q  (kW/m2) 

481.4 

477.7 

409.6 

h0  (kW/m2K) 

9.86 

9.82 

10.75 

AT(K) 

48.82 

48.65 

38.10 

Figure  8  shows  the  wall  temperature  profiles  for  the  data  points  listed  in  Table 
3  along  with  a  vacuum  data  set  (P^—28  kPa,  Vvapor~l  m/s,  V^^—2.75  m/s)  for 
comparison,  where  0  degrees  is  at  the  top  dead  center  position  of  the  tube.  The 
shape  of  the  temperature  profiles  shows  the  effect  of  condensate  film  thickening 
toward  the  bottom  of  the  tube.  The  higher  resistance  through  a  thicker  condensate 
film  results  in  a  lower  tube  wall  temperature  toward  the  bottom  of  the  tube  as 
shown. 

As  expected,  the  wire  wrap  insert  and  Heatex  insert  temperature  profiles  at 
atmospheric  pressure  are  very  similar.  The  no  insert  profile  has  a  shape  similar  to 
the  wire  wrap  and  Heatex  profiles,  but  the  mean  wall  temperature  is  about  UK 
higher  and  shows  the  effect  of  the  increased  inside  resistance  to  heat  transfer.  The 
vacuum  run  temperature  profile  shows  the  effect  of  a  lower  temperature  gradient 
between  the  steam  and  coolant;  the  lower  heat-transfer  potential  results  in  lower 
heat  fluxes  and  a  flatter  temperature  profile. 

Figure  9  shows  the  comparison  of  the  instrumented  data  with  the  predictions 
of  Nusselt,  Fujii,  and  Shekriladze  and  Gomelauri  covered  in  section  II  B.  The  data 
depicted  ranges  from  Psal  =  28  kPa  and  Vvapor  =  3.5  m/s  to  Psat  =  101  kPa  and  Vvapor 
=  1.1  m/s.  The  data  seems  to  follow  the  Shekriladze  and  Gomelauri  prediction  the 
closest,  but  is  also  very  near  the  Fujii  prediction. 
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Figure  8.   Horizontal  Tube  Wall  Temperature  Profiles  (Smooth 
Instrumented  Tube,  Vvapor~l  m/s) 
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Figure  9.   Comparison  of  Experimental  Results  with  the  Predictions 
of  Nusselt,  Fujii,  and  Shekriladze-Gomelauri. 
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C.       INSIDE   HEAT-TRANSFER   CORRELATIONS    FROM    INSTRUMENTED 
TUBE  RESULTS 

The  instrumented  tube  provides  a  method  of  determining  an  inside  heat- 
transfer  coefficient  for  each  insert  from  direct  wall  temperature  measurements. 
Wanniarachchi  et  al  [Ref.  8],  in  a  similar  effort  to  resolve  the  differences  between 
the  Sieder-Tate  correlation  and  his  experimental  results,  developed  a  correlation 
based  on  the  Sieder-Tate  expression.  Using  an  intercept  form,  with  data  taken  on 
the  same  apparatus  as  used  in  this  study,  and  a  least-squares  fitted  leading 
coefficient,  the  correlation  took  the  following  form: 


Nu  =0.064  ReoiPrm 


-^       +26.4  (6.1) 


However,  Rouk  [Ref.  4],  using  an  optimization  technique,  showed  that  the  value  of 
the  intercept  had  little  effect  on  the  results,  and  it  was  the  accurate  determination 
of  the  Reynolds  exponent  that  was  more  critical. 

With  a  view  to  finding  the  appropriate  Reynolds  exponent  and  leading 
coefficient  for  equation  (5.23),  the  data  for  no  insert,  wire  wrap,  and  Heatex  insert 
were  plotted  as  In  Re  versus  In  (NufPrm (jijfiw)° 14)  as  explained  in  section  V  C.  The 
plotted  data  are  shown  in  Figures  10,  11  and  12  for  no  insert,  wire  wrap  insert,  and 
Heatex  insert  respectively.  A  line  of  best  fit  (typically  with  a  regression  coefficient 
of  0.99)  was  used  to  obtain  the  value  of  the  intercept.  Figure  13  shows  the  plot  of 
all  three  cases  on  the  same  graph  for  comparison.  Note  that  the  insert  data  is 
closely  grouped  on  the  upper  regions  of  the  graph  when  compared  to  the  no  insert 
data.  The  increase  in  the  Nusselt  number  again  indicates  more  efficient  inside  heat 
transfer  for  the  insert  vice  the  no  insert  case.    The  Reynolds  exponent,  or  slope, 


57 


m 


Figure  10.   Log-Log  Plot  of  Re  versus  Nu/Pr173^,//^)014  for  No 
Insert  (Smooth  Instrumented  Tube) 


58 


m 


Figure  11.   Log-Log  Plot  of  Re  versus  Nu/Pr1/3(^«//^)° 14  for  Wire  Wrap 
Insert  (Smooth  Instrumented  Tube) 
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Figure  12.   Log-Log  Plot  of  Re  versus  Nii/Pr^O^)014  for  Heatex 
Insert  (Smooth  Instrumented  Tube) 
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Figure  13.   Combined  Log-Log  Plot  of  Re  versus  Nuf?Tm(aJn^) 
Insert  Conditions  (Smooth  Instrumented  Tube) 
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differs  from  the  Sieder-Tate-type  equation  of  0.8  in  all  three  cases.  The  following 
derived  correlations  apply  specifically  only  to  the  medium  tube,  but  should  be 
applicable  to  any  tube  with  the  same  inside  diameter. 

The  no  insert  inside  heat-transfer  correlation  had  the  form: 


Nu  =0.013  Re°*9Prir* 


I  u    \°  14 


^c 


\Pw/ 


(6.2) 


Equation  (6.2)  was  used  to  reprocess  some  current  smooth  tube  runs,  with  no  insert, 
as  well  as  those  of  Guttendorf  [Ref.  3]  and  Van  Petten  [Ref.  2]  to  check  their  values 
for  the  smooth  tube  Nusselt  coefficient  a. 

The  wire  wrap  insert  inside  heat-transfer  correlation  had  the  form: 


/  ii  \014 


Nu  =0.052  Re°*2Priri 


(6.3) 


Equation  (6.3)  was  used  to  reprocess  previous  data  taken  on  the  same  apparatus  by 
Guttendorf  and  Van  Petten,  who  used  the  wire  wrap  insert  for  their  finned  and 
smooth  tube  experiments,  with  a  view  to  checking  their  results  with  this  new 
correlation.  Table  4  shows  a  comparison  between  the  new  wire  wrap  insert 
correlation  (eq.  6.3),  and  Wanniarachchi's  correlation  (eq.  6.1)  both  with  and  without 
the  intercept  value  included.  This  comparison  shows  that  equation  6.1  is  actually 
more  accurate  without  the  intercept  value  included  when  compared  to  the  results  of 
the  new  wire  wrap  insert  correlation. 
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Table  4.   COMPARISON  OF  EQUATIONS  (6.1)  AND  (6.3) 

FOR  Pr1/3(/i<//O014=l-4  (held  constant  for  comparison) 

Reynolds  Number 

20,000 

30,000 

40,000 

50,000 

Nusselt  Number 
Nu=0.064  Re08Pr1/3(/x,//O014+26.4 

273.7 

368.4 

456.9 

541.0 

Nusselt  Number 
Nu=0.064  Re08Pr1/3(/axw)014 

247.3 

342.0 

430.5 

514.6 

Nusselt  Number 
Nu= 0.052  Re°-82Pr1/3(/A JiiJu 

244.9 

341.5 

432.3 

519.1 

The  Heatex  insert  inside  heat- transfer  correlation  had  the  form: 


Nu=022  Re069Pr113 


(*c 


\0.14 


I,  ^"/ 


(6.4) 


Equation  (6.4)  was  used  to  reprocess  all  Heatex  data. 

Memory  [Ref.  12]  conducted  condensation  experiments  on  a  different  apparatus 
with  a  smooth  instrumented  tube.  He  also  determined  the  Reynolds  number 
exponent  for  no  insert,  wire  wrap  insert  (made  locally  and  somewhat  different  from 
the  one  used  in  this  study),  and  Heatex  insert.  The  exponents  he  obtained  are 
reported  in  Table  5.  The  Heatex  insert  exponent  of  0.68  compares  very  favorably 
with  the  value  of  0.69  found  in  this  study.  The  wire  wrap  exponent  of  0.73  was  well 
below  the  value  of  0.82,  most  likely  due  to  differences  in  the  insert.  The  no  insert 
case  gave  a  value  of  0.85  and  compared  well  to  the  value  of  0.89  found  in  this  study. 
The  difference  in  the  no  insert  exponent  is  most  likely  due  to  the  difference  in  the 
tube  entrance  region;  Memory  used  a  long  run  of  straight  pipe  for  the  tube  entrance, 
whereas  the  present  work  had  a  sharp  bend  just  prior  to  the  condenser  tube. 

Rouk  [Ref.  4]  used  the  instrumented  smooth  tube  data  of  Georgiadis  [Ref.  5] 
to  the  find  the  appropriate  value  of  the  Reynolds  number  exponent  for  the  no  insert 
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and  wire  wrap  cases.    His  results  compare  quite  well  with  this  study  and  are  also 
given  in  Table  5. 

ANL  (Argonne  National  Laboratory)  [Ref.  21]  conducted  an  assessment  of 
heat-transfer  correlations  for  turbulent  pipe  flow  with  water  to  determine  the  best 
correlation (s)  on  which  to  base  their  design  of  Ocean  Thermal  Energy  Conversion 
(OTEC)  heat  exchangers.  ANL  used  two  shell-and-tube  heat  exchangers,  with  no 
inserts,  for  analysis  and  reported  the  following: 

1.  The  Dittus-Boelter  (eq  2.7)  and  Sieder-Tate  (eq  2.9)  correlations  under- 
predicted  the  data  by  5%  to  15%  and  were  considered  too  conservative  for 
design. 

2.  Overall,  the  "best"  correlations  were  found  to  be  Petukhov-Popov  (eq  6.5)  and 
Sleicher-Rouse  (eq  6.6),  both  of  which  showed  excellent  agreement  (±  5%) 
with  the  experimental  data  (at  Pr=6.0  and  Pr  =  11.6). 


Nu  = (e,S)RePr 

#1  +  £2(e/8)1/2(/V2/3-l) 

(valid  for  0.5  <  Pr  <  2000  and  104  <  Re  <  5xl06) 


where: 

€ 

K,      =    1  +  3.4€ 


€        =    (1.82  log10  Re  -  1.64) 2 


-1/3 


K2      =    11.7  +  1.8  Pr 

Nu  =  5  +0.015  RefPr*  (6-6) 

(valid  for  0.1  <  Pr  <  105  and  104  <  Re  <  106) 


where: 


a      =    0.88  -  0.24/(4+Prw) 
b      =    1/3  +  0.5e*6Plw 
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3.  The  most  accurate  correlations  (i.e.  Petukhov-Popov  and  Sleicher-Rouse) 
seem  to  yield  effective  Reynolds  exponents  in  the  neighborhood  of  0.85 
(uncertainty  range:  m  =  0.82  to  0.88). 

4.  The  potential  sources  of  uncertainty  in  the  Wilson  procedure  included 
waterside  flow  maldistribution,  entrance  effects,  experimental  error  in  U  ,  and 
the  uncertainty  in  the  Reynolds  number  exponent.  Of  these,  they  concluded 
that  the  uncertainty  in  the  Reynolds  number  exponent  was,  by  far,  the  most 
significant.  In  fact,  the  results  of  the  Wilson  procedure  were  found  to  be 
highly  sensitive  to  the  value  of  the  Reynolds  number  exponent. 


Table  5. 

COMPARISON  OF  REYNOLDS  NUMBER  EXPONENTS  FOR 
SIEDER-TATE-TYPE  CORRELATIONS 

Experimental 
Data 

Rouk 

Memory 

ANL 

No  Insert 

0.89 

0.90 

0.85 

0.85 

Wire  Wrap 

0.82 

0.78 

0.73 

— 

Heatex 

0.69 

— 

0.68 

— 

D.      ANALYSIS  OF  SMOOTH  TUBE  RESULTS 

When  using  the  modified  Wilson  plot  technique  to  reprocess  data  files,  the 
solution  option  can  be  specified  to  use  either  the  stored  value  of  the  Sieder-Tate 
coefficient  (for  direct  computation  of  h0)  or  let  the  coefficient  value  "float",  which 
allows  the  program  to  calculate  its  own  value  of  the  coefficient.  In  order  to 
determine  which  method  was  most  accurate,  the  instrumented  data  files  were 
reprocessed  using  each  method  and  then  compared  with  the  values  of  the  heat- 
transfer  coefficient  which  were  obtained  by  direct  measurement  of  the  tube  wall 
temperature.  A  high,  medium,  and  low  coolant  flow  rate  was  chosen  from  each  run 
to  facilitate  the  comparison.  The  results  were  tabulated  and  are  shown  in  Appendix 
D;  it  can  be  seen  that  the  fixed  coefficient  method  yielded  the  more  accurate  results 
at  least  75%  of  the  time.  The  mean  error  of  the  fixed  method  was  ±2.0%,  and  that 
of  the  floating  method  was  ±5.4%.   The  error  for  the  lowest  coolant  flow  rate  (Re 
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<  20,000)  was  noticeably  higher  than  for  higher  coolant  flow  rates  for  both  methods. 
The  choice  of  using  the  fixed  coefficient  method  represents  a  departure  from  the 
practice  of  previous  researchers  on  this  apparatus  who  exclusively  used  the  floating 
coefficient  method. 

Prior  to  the  instrumented  tube  runs  reported  in  Figure  5  (from  which  the  new 
correlations  were  empirically  derived)  a  series  of  data  runs  were  made  using  a  plain 
smooth  uninstrumented  tube  (S02)  of  the  same  dimensions  using  the  Heatex  insert. 
The  plain  smooth  tube  data  was  then  reprocessed  using  the  new  correlation  for  the 
Heatex  insert  and  the  results  are  plotted  in  Figure  14.  These  data  sets  were  taken 
at  the  same  conditions  as  the  instrumented  data  of  Figure  5  except  for  the  set  at  the 
highest  vapor  velocity  of  6.2  m/s  vice  3.5  m/s  for  the  instrumented  tube.  Similar 
effects  of  vapor  shear  and  vapor  pressure,  as  mentioned  previously  for  Figure  5,  are 
clearly  seen,  and  again  illustrate  the  vapor  shear  effect  on  the  outside  heat-transfer 
coefficient. 

With  the  exception  of  the  two  data  runs  at  high  vapor  velocity,  the  data  from 
Figures  5  and  14  are  shown  together  in  Figure  15.  The  close  agreement  of  the 
reprocessed  plain  smooth  tube  data  with  the  instrumented  smooth  tube  data  allows 
a  high  degree  of  confidence  in  the  accuracy  of  the  new  correlations  and  the  choice 
of  the  fixed  coefficient  method. 

To  provide  a  baseline  from  which  to  evaluate  finned  tube  performance  it  was 
necessary  to  obtain  the  smooth  tube  Nusselt  coefficient,  a,  for  the  specific  conditions 
under  which  the  comparison  was  to  be  made.  The  condition  chosen  was  atmospheric 
pressure  (101  kPa)  and  a  vapor  velocity  of  ~1  m/s. 

For  8  complete  sets  of  data  the  average  value  of  a,  using  the  fixed  method,  was 
found  to  be  0.876.   The  average  value  for  each  data  set  was  found  by  taking  the 
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Figure  14.   Effect  of  Pressure  and  Vapor  Velocity  on  the  Steam 
Heat-Transfer  Coefficient  (Non-instrumented  Smooth 
Tube,  Heatex  Insert) 
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Figure  15.   Comparison  of  Instrumented  Smooth  Tube  Results  with  Non- 
Instrumented  Smooth  Tube  Data  After  Reprocessing  with  the 
New  Heatex  Insert  Inside  Heat-Transfer  Correlation. 
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measured  value  of  h0  for  each  data  point,  dividing  it  by  the  Nusselt  theory  prediction 
of  h0,  and  then  multiplying  by  0.728  (the  Nusselt  coefficient);  the  average  of  all  the 
data  points  in  the  set  was  then  taken.  Interestingly  the  Wilson  plot  floating 
coefficient  method  gave  a  value  of  0.835,  somewhat  lower  than  the  average  value. 
Originally  it  was  thought  that  this  discrepancy  might  be  due  to  "outlier"  data  points 
(high  or  low  coolant  velocity)  in  each  set.  However,  removing  the  highest  or  lowest 
coolant  flow  rates  within  a  set  had  little  or  no  effect  on  the  Wilson  plot  result.  The 
reason  for  the  discrepancy  is  still  not  known  and  merits  future  study. 

The  value  of  a  was  calculated  for  several  other  flow  conditions;  these  are 
shown  in  Table  6.  The  trend  of  the  readings,  like  that  of  Figure  6,  shows  that  vapor 
velocity  has  a  much  greater  effect  than  pressure  on  the  value  of  a,  as  expected. 

Table  6.  SMOOTH  TUBE  a  SUMMARY;  EFFECT  OF  PRESSURE  AND 
VAPOR  VELOCITY 


File  Name 

P(kPa) 

Vapor  Velocity 

(m/s) 

a 

FIMAVSH1 

28 

3.5 

1.015 

FIMAVSH2 

41 

2.5 

.985 

FIMAVSH3 

68 

1.5 

.930 

FIMAVSH4 

101 

1.1 

.866 

FIMAVSH7 

69 

1.0 

.836 

FIMAVSH6 

41 

1.0 

.818 

FIMAVSH5 

28 

1.0 

.786 

E.      ANALYSIS  OF  FINNED  TUBE  RESULTS 

With  an  accurate  value  of  a,  and  the  newly  determined  inside  heat-transfer 
correlations,  the  medium  family  finned  tube  data  of  Van  Petten  [Ref.  2]  was 
evaluated.  Figure  16  shows  the  data  Van  Petten  reported  in  his  thesis;  it  also  shows 
his  data  after  being  reprocessed  using  the  new  wire  wrap  insert  correlation  with 
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Figure  16.   Comparison  of  the  Steam  Heat-Transfer  Enhancement  Data  of  Van 
Petten,  for  the  Medium  Finned  Tube  Family,  Using  the  Modified 
Wilson  Plot  and  New  Wire  Wrap  Insert  Inside  Correlation. 


70 


both  the  fixed  coefficient  and  Wilson  plot  floating  coefficient  methods.  Since  Van 
Petten  used  the  Wilson  plot  method,  it  is  not  surprising  that  the  original  thesis  data 
and  new  Wilson  plot  floating  coefficient  data  are  comparable  since  the  Reynolds 
exponent  only  varied  from  0.8  to  0.82.  The  fixed  coefficient  method  enhancement 
is  substantially  higher  than  the  Wilson  plot  results. 

Since  the  assertion  is  that  the  fixed  coefficient  method  is  more  accurate  than 
the  Wilson  plot  method,  then  the  conclusion  must  be  that  the  enhancement  for  this 
set  of  finned  tubes  is  actually  higher  than  previously  reported. 

During  this  study,  limited  medium  finned  tube  experiments  were  conducted  for 
purposes  of  comparison.  Figure  17  shows  the  comparison  between  this  data  and  the 
newly  reprocessed  data  of  Van  Petten  (using  the  fixed  coefficient  method)  and  shows 
reasonable  agreement.  To  more  clearly  illustrate  this  point,  the  data  taken  on  the 
2.0  mm  fin  spacing  tube  has  been  given  in  more  detail  in  Figure  18.  Excellent 
agreement  is  seen  between  the  experimental  results  of  this  study  and  that  of  Van 
Petten  using  the  known  inside  heat-transfer  correlation  with  the  fixed  coefficient 
method.  Again,  as  shown  in  Figure  16,  the  Wilson  plot  prediction  is  significantly 
below  the  fixed  coefficient  results. 
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Figure  17.   Comparison  of  the  Steam  Heat-Transfer  Enhancement  Data 

of  Van  Petten  and  Swensen  for  the  Medium  Finned  Tube  Family. 
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Figure  18.   Comparison  of  the  Steam  Heat  Transfer  Data  of  Van  Petten  and 
Swensen  for  the  2.0  mm  Fin  Spacing  Medium  Tube. 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 
A.       CONCLUSIONS 


1.  The  inside  heat- transfer  correlation  is  highly  sensitive  to  the  Reynolds  number 
exponent. 

2.  Each  insert  condition  must  be  analyzed  separately  to  determine  the 
appropriate  "form"  of  the  inside  heat- transfer  correlation. 

3.  Calculations  based  on  a  known  inside  heat-transfer  correlation  are  more 
accurate  than  modified  Wilson  plot  results. 

4.  Armed  with  accurate  inside  heat-transfer  correlations,  previous  data  may  be 
reprocessed  to  give  more  accurate  results. 

5.  The  source  of  contamination  in  the  test  apparatus,  which  has  caused  a 
dropwise  condensation  problem  for  a  number  of  years,  is  most  probably  due 
to  a  contaminated  distilled  water  source. 


B.       RECOMMENDATIONS 


1.  Reprocess  all  previous  medium  and  large  diameter  finned  tube  data  using  the 
fixed  coefficient  method  to  obtain  more  accurate  results. 

2.  Continue  with  construction  of  smooth  instrumented  tubes  of  different 
diameters  (i.e.  small,  medium,  and  large)  to  confirm  the  medium  tube  results 
and  develop  correlations  specifically  for  the  small  and  large  diameter  tubes. 

3.  Construct  one  representative  instrumented  finned  tube  to  test  the  validity  of 
applying  instrumented  smooth  tube  results  to  finned  tube  data. 

4.  Test  representative  water  samples  that  have  been  collected  from  both  the 
distiller  and  boiler  to  confirm  the  presence  of  impurities  and  validate  their 
origin. 

5.  Replace  current  distiller  with  a  deionized  pure  water  source  (either 
commercial  purchase  of  water  or  new  distilling  apparatus). 
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APPENDIX  A.    PHYSICAL  AND  THERMODYNAMIC  PROPERTIES  OF 
WATER 


The  physical  and  thermodynamic  properties  of  water  were  based  on  the 
following  equations: 

Heat  Capacity  (J/kg-K): 

Cp  =  4.211  -r*[2.268xl0-3-r*(4.424^10-5+2.714xl0-7*  7)]  (A-l) 

where:  T  =  temperature  (celsius) 

Dynamic  viscosity,  (kg/ms): 

u.  =  (2 Ax  10"5)  *  ltf247  W 13315>]  (A.2) 

Thermal  conductivity  (W/m-K): 

k  =  -0.9225+x  *(2.8395  -x  *(1.8007  -x  *(0.5258  -0.0734  **)))  (A3) 

where:  x  =  (T+273.15)/273.15 

Density  (kg/m3): 
p  =999.5295  +  T*  (0.0127 -T*  (5.4825x70  3-T*  1.2341x70  5))  (A-4) 
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Latent  heat  of  vaporization  (J/kg): 

hh  =  2477200  -  2450  *  (T- 10)  ( A.5) 

Fluid  Enthalpy  (J/kg): 
hf=  T*  (4.2038  -  T*  (5.8813x70-4-r*4.5516jc70-6))  (A.6) 


Saturation  pressure  (Pa): 


P    =22120000  *Pr 

sat 

where: 

Pr 

=    eBr 

Br 

=  suM/n 

fr(l+  4.1671  (1-T)  +  20. 

(A.7) 


=    SUM  /[Tr  (1+  4.1671  (1-T)  +  20.9751  (1-Tr)2)  -  (1-Tr)/ 
((1.0xL10Wr)2  +  6V)] 

SUM    =  (-7.6912)(1-Tr)  -  26.0802  (1-Tr)2  - 168.1707  (1-Tr)3  +  64.2329  (1-Tr)4 
118.9646  (1-Tr)5 

Tr         =  (T+273.15)/  647.3 
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APPENDIX  B.    SYSTEM  CALIBRATIONS  AND  CORRECTIONS 
B.l    Thermocouple  and  Quartz  Thermometer  Calibration 

Several  different  thermocouples  and  the  quartz  crystal  thermometer  (HP 
2804 A,  Ser.  No.  2244AD1192)  were  calibrated  against  a  platinum  resistance  probe 
in  a  mixed  isothermal  ethylene  glycol  bath  with  a  Rosemont  Galvanometer  model 
920A  commutating  bridge  (Ser.  no.  013494)  from  23  to  26  September,  1991. 

NOMENCLATURE 

Tl  Quartz  Thermometer  Measuring  Probe  Tl  (Ser.  No.  2120A-00707,  with  dial 

setting  481) 
T2  Quartz  Thermometer  Measuring  Probe  T2  (Ser.  No.  2120A-60459  with  dial 

setting  510) 
T1-T2         The  T1-T2  reading  on  the  quartz  thermometer 
T-55  Large  diameter  metal  sheath  thermocouple  (diameter  =  0.040") 

T-56  Small  diameter  metal  sheath  thermocouple  (diameter  =  0.020") 

T-57  Old  Thermocouple  (taken  from  rig  during  disassembly  set  to  HP  3497 A 

internal  zero);  (Type:  Omega,  TT-T-30,  Lot#  HCP093HC0306) 
T-58  New  1  Thermocouple  (set  to  HP  3497 A  internal  zero);  (Type:  Omega,  TT- 

T-30  SLE,  Lot#  OCP1453PTCC01473P) 
T-59  New  2  Thermocouple  (referenced  to  ice  bath  zero) 

T-60  10  Junction  Thermopile  #1  (referenced  to  ice  bath  zero) 

T-61  10  Junction  Thermopile  #2  (referenced  to  ice  bath  zero) 

T-62  10  Junction  Thermopile  #3  (referenced  to  ice  bath  zero) 

The  calibration  was  performed  by  taking  the  bath  temperature  up  from  290  K 
to  393  K  by  5  K  increments  then  back  down  to  check  for  hysteresis;  no  hysteresis  was 
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observed.  The  data  was  fitted  to  a  fifth-order  polynomial  (regression  coefficient  = 
1.000  for  each  polynomial  fit)  in  each  case  with  the  following  results: 

T55  =  273.15+(2.5943e-2)K-(7.2671-7)K2 

+  (3.2941-11)F3-(9.7119-16)K4  (B.l) 

+  (9.7121 -20)** 


T56  =  273.15+(2.5878€-2)K-(5.9853e-7)K2 

-  (3.1242£-ll)K3+(1.3275e-14)K4  (B-2) 

-  (1.0188^-18)^ 


757  =  273.15+(2.5923e-2)K-(7.3933e-7)K2 

+  (2.8625*-ll)73+(1.97i7e-15)K4  (B3) 

-  (2.2486e-19)K5 


T58  =  273.15+(2.5931*-2)K-(7.5323e-7)K2 

+  (4.0567e-ll)K3-(1.2791e-15)K4  (B4) 

+  (6.4402^-20)^ 


759  =  273.15+(2.5471e-2)K-(3.7621e-7)K2 

-  (1.0105e-10)K3+(2.3928e-14)K4  (B5) 

-  (1.6440^-10)^ 


T60  =  273.17+(2.5571e-2)K-(1.9980e-7)K2 

-  (1.6385*-10)K3+(2.6164€-14)K4  (B6) 

-  (1.0295^-18)^ 


T61  =  273.15+(2.6119e-2)K-(9.0449*-7)K2 

+  (1.1214*-10)K3-(1.5623<?-14)F4  (BJ) 

+  (1.0646^-18)^ 
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7152  =  273.15+(2.5996«-2)K-(7.4405e-7)K2 

+  (2.4733e-ll)K3+(3.323fc?-15)K4  (B-8) 

-  (3.7460e-19)K5 

For  T55,  T56,  T57,  T59;  V=  Voltage  in  microvolts  (eg.  0.010  volts  =  10,000 
microvolts). 

For  T60,  T61,  T62;  V=  Voltage  in  microvolts/10  (eg.  100,000  microvolts  +  10 
=  10,000). 

The  HP  2804A  Quartz  Crystal  Thermometer  was  also  calibrated  on  26 
September,  1991  with  the  results  summarized  in  Table  B.l. 


Table  B.l. 

QUARTZ  THERMOMETER  CALIBRATION  DATA 

Reference 

Temperature 

(deg  C) 

77 
(deg  C) 

Error 
(deg  C) 

T2 
(deg  C) 

Error 
(decC) 

17.979 

18.003 

0.024 

18.001 

0.022 

18.977 

18.997 

0.020 

18.996 

0.021 

19.984 

20.001 

0.017 

20.002 

0.018 

20.980 

21.007 

0.027 

21.008 

0.028 

22.936 

22.959 

0.023 

22.961 

0.025 

25.014 

25.032 

0.018 

25.037 

0.023 

30.140 

30.153 

0.013 

30.164 

0.024 

34.977 

34.983 

0.006 

34.998 

0.021 

The  mean  error  for  the  Quartz  Thermometer  as  shown  in  the  Table  B.l  is 
— +0.02  deg  C.  However,  it  is  well  noted  that  T1-T2,  the  critical  measurement, 
indicates  an  apparent  error  of  less  than  +0.005  deg  C  when  the  temperature  reading 
is  near  or  below  25  deg  C;  77  and  72  tracking  well  together  lowers  the  error 
estimate. 
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B.2   Flow  Meter  Calibration 

The  flow  meter  calibration  for  coolant  flow  through  the  single  horizontal  tube 
was  completed  on  Oct  28,  1991  using  a  stop  watch,  portable  tank,  and  a  Toledo 
model  31-0851  IV,  Se.  No.  1326  scale  with  1/10  pound  graduations.  The  following 
relation  was  obtained  via  linear  regression: 


m  =  6.7409F  +  13.027 

B.9 


where: 
m 


=    mass  flow  rate  (grams  per  second) 
=    flow  meter  reading  (eg.  10%=>10) 


The  applicable  range  of  the  calibration  was  10%  to  95%.  The  water 
temperature  on  the  day  of  the  calibration  was  17.5°C  (290.6  K)  and  water  density  was 
998.5  kg/m3.   The  data  is  shown  in  Figure  B.l 

B.3   Mass  Flow  Rate  Correction 

The  inlet  water  temperature  from  the  cooling  water  sump  varies  anywhere  from 
15°C  to  25°C  depending  on  environmental  conditions.  To  account  for  these 
temperature  variations  the  following  function  was  used  to  calculate  a  correction 
factor  for  viscosity  variation  with  temperature  [Ref.  5]. 

C/=1.0365  -(1.9644£-3)77n  +(5.2500e-6)7w2  (B1°) 


where: 


Cf     =    mass  flow  rate  correction  factor 
Tin     =    inlet  temperature  (celsius) 
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Figure  B.l   Horizontal  Tube  Coolant  Flowmeter  Calibration  Chart 
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The  value  of  Cf  for  the  flow  meter  calibration  (Tw  =  17.5°C)  was  1.0037. 
Therefore,  the  actual  mass  flow  rate  was  calculated  using  the  following  equation: 

m     =m  .-£-  (B.ll) 

"*      "*  1.0037 

where: 

mact    =    corrected  mass  flow  rate  (grams  per  second) 

mca,c=    computed  mass  flow  rate  (eq.  B.9)  (grams  per  second) 

B.4    Pressure  Transducer  Calibration 

Three  methods  of  pressure  measurement  were  available  on  the  apparatus: 

1.  Direct  pressure  reading  off  the  Heise  solid  front  -  CM-1041 19  pressure  gauge, 
(range  0-15  psia). 

2.  Converted  voltage  readings  from  the  Setra,  model  204,  Ser.  no.  63982 
pressure  transducer  (range  0-14.7  psia;  0-5  volts;  5V— 0  psia). 

3.  Steam  saturation  temperature  measurement  with  the  apparatus  producing 
steam  at  steady  state.  The  steam  saturation  temperature/pressure  relation 
was  utilized  via  standard  steam  tables. 

The  pressure  transducer  was  calibrated  versus  the  vapor  temperature  probe 

reading  on  12  December  1991.   Equation  B.12  gives  the  desired  relationship.   The 

data  is  shown  in  Figure  B.2. 

P=-2.9360K+ 14.7827  (B12) 


where: 


P        =    pressure  (psia) 

V       =    pressure  transducer  voltage  reading  (volts) 


B.5    FRICTION  TEMPERATURE  CORRECTION 

As  coolant  flows  through  the  tube  there  is  a  bulk  temperature  rise  due  to 
frictional  heating,  which  is  highly  dependent  on  fluid  velocity.   Although  small,  this 
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Figure  B.2   Pressure  Transducer  Calibration  Chart 
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can  have  a  significant  effect  on  the  calculated  overall  heat-transfer  coefficient. 
Measurements  were  made  for  no  insert,  Heatex,  wire  wrap  and  twisted  tape  inserts 
as  shown  in  Figure  B.3,  on  5  December  1991.  Each  data  set  was  curve  fitted  to  a 
third  order  polynomial  which  is  depicted  in  Table  B.2.  Each  respective  polynomial 
corrects  the  temperature  rise  measurement  for  the  heating  due  to  the  particular  type 
of  insert  used. 

Table  B.2   FRICTION  TEMPERATURE  RISE  POLYNOMIALS 


Insert  Type 

Polynomial 

None 

Trise  =  -1.960x10  5V3+9.349xl0"4V2+  1.749x10  4V-2.728xl0"4 

Wire  Wrap 

Trise  =  8.160x10  5V3+ 1.4512x10  3V2+2.745xl03V-3.991xl0"4 

Heatex 

Trise  =  8.160x10  5V3+ 1.080x10  3V2+ 1.232x10  3V+ 8.570x10 5 

Twisted  Tape 

Trise  =  4.070x10  5V3+4.451xl0"4V2+ 1.711x10  3V-6.440xl05 

where:  Trise  =  temperature  rise  (K) 
V  =  fluid  velocity  (m/s) 

84 


tr 

0) 

CO 

tr 

0 
(0 

CD 

c 

Q. 

0) 

— 

cn 

CO 

Q. 

tr 

I— 

c 

CT3 

CD 
CO 

CD 

X 

£ 

c 

CD 

CO 

CT? 

CD 

o 

1- 

CD 

IE 

£ 

CD 


-  in 


-  ^ 


-  CO 


(0 

E 


o 
o 

> 

o 


-  CM 


fo 


00 
O 

o 

O 

lO 
O 

O 

CO 

o 

CM 
O 

O 

o 
o 

O 

o 

O 

O 

O 

o 

O 

O 

o 

(X)    asm    ejnjBJeduiei 


Figure  B.3   Friction  Temperature  Rise  Curves  for  Heatex  Insert,  Wire  Wrap 
Insert,  Twisted  Tape  Insert,  and  No  Insert. 
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APPENDIX  C.    SYSTEM  INTEGRITY  /  LEAK  TESTING 

As  mentioned  in  section  III  C,  the  material  for  the  auxiliary  condenser 
penetration  plates  was  changed  from  aluminum  to  stainless  steel  on  24  January  1992. 
The  stainless  steel  screw  thread  connectors  in  the  aluminum  side  plate  had  loosened 
to  the  extent  that  leakage  could  be  detected.  The  cause  of  this  loosening  was  due 
to  thermal  cycling  of  the  apparatus  and  the  differential  contraction/expansion  of  the 
aluminum/stainless  steel  combination. 

An  initial  leak  test  was  conducted  from  20  December  1991  through  2  January 
1992;  the  results  are  shown  in  Figure  C.l.  The  initial  mean  leak  rate  was  3.4  mmHg 
per  day. 

Subsequent  to  the  structural  modification  noted  above,  another  leak  test  was 
conducted  6-19  February  1992;  the  results  are  shown  in  Figure  C.l.  The  mean  leak 
rate  was  found  to  be  1.7  mmHg  per  day,  an  noticeable  improvement.  In  general,  a 
leak  rate  of  2  mmHg  per  day  is  considered  acceptable. 
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Figure  C.l   Apparatus  Leak  Test  I 
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Figure  C.2   Apparatus  Leak  Test  II 
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APPENDIX  D.    COMPARISON  OF  FIXED  C,  vs  FLOATING  C,  SOLUTION 
METHODS  FOR  MODIFIED  WILSON  PLOT  DATA  REPROCESSING. 


As  related  in  section  VI  D  the  "fixed"  C,  method  and  "floating"  C,  method  were 
evaluated  against  the  original  instrumented  tube  results  in  order  to  choose  the  most 
accurate  method  of  reprocessing  non-instrumented  data.  The  results  of  the  fixed 
coefficient  versus  modified  Wilson  plot  floating  coefficient  method  comparison  are 
shown  in  Table  D.  The  fixed  coefficient  method  was  determined  to  be  more  accurate 
than  the  floating  coefficient  method  75%  of  the  time.  For  the  three  coolant  flow 
rates  considered  for  each  run,  the  overall  fixed  coefficient  method  mean  error  was 
±2.0%,  and  the  floating  coefficient  method  mean  error  was  ±5.4%. 

For  Re  < 20,000  the  error  was  noticeably  higher,  which  tends  to  support  the 
assertion  of  a  number  of  researchers  that  data  with  Reynolds  numbers  below  20,000 
should  not  be  used,  since  the  flow  may  not  be  fully  turbulent. 
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Table  D.  COMPARISON  OF  FIXED  C,  vs  FLOATING  C,  REPROCESSING 
METHODS 


File  Name 

Flow 
Vel. 

Film 
AT 

Inst 
Tube 
Data 

Fixed 
Method 

% 
error 

Float 
Method 

% 
error 

vw 

(m/s) 

k 

h°2 

(kW/m2K) 

h0 

h0 

FIMAVSH1 

4.35 

29.33 

12.13 

11.93 

-1.7 

13.19 

8.7 

2.21 

22.81 

13.18 

12.82 

-2.6 

15.17 

15.1 

1.16 

18.27 

13.68 

13.16 

-3.8 

-2.7* 

17.20 

25.7 
+  16.5 

FIMAVSH2 

4.34 

36.09 

11.10 

11.14 

0.4 

11.07 

0.3 

2.22 

30.39 

12.03 

11.87 

-1.3 

11.72 

-2.6 

1.16 

24.62 

12.66 

12.50 

-1.3 

12.26 

-3.1 

-0.7*  -1.8 


FIMAVSH3 

4.35 

46.10 

10.31 

10.26 

-0.5 

10.62 

3.0 

2.23 

39.35 

10.89 

10.83 

-0.6 

11.44 

5.1 

1.16 

32.56 

11.42 

11.40 

-0.2 
-0.4* 

12.43 

8,8 
+5.6 

FIMAVSH4 

4.35 

55.20 

9.263 

9.146 

-1.3 

9.584 

3.5 

2.23 

47.31 

9.869 

9.756 

-1.1 

10.52 

6.6 

1.16 

39.45 

10.41 

10.39 

-0.2 
-0.9* 

11.71 

12.5 

+7.5 

FIMAVSH5 

4.35 

32.83 

8.618 

8.478 

-1.6 

9.018 

4.6 

2.23 

28.05 

9.589 

9.341 

-2.6 

10.36 

8.0 

1.16 

22.99 

10.23 

9.865 

-3.6 

11.66 

14.0 

-2.6*  +8.9 
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File  Name 

Flow 
Vel. 

Film 
AT 

Inst 
Tube 
Data 

Fixed 
Method 

% 
error 

Float 
Method 

% 
error 

vw 

(m/s) 

% 

h°2 
(kW/m2K) 

ho 

ho 

FIMAVSH6 

4.35 

39.23 

9.073 

8.953 

-1.3 

9.084 

0.1 

2.23 

33.76 

9.757 

9.447 

-1.3 

9.667 

1.0 

1.16 

27.98 

10.18 

10.25 

06 
-0.7* 

10.64 

15 
+  1.9 

FIMAVSH7 

4.35 

49.04 

8.960 

8.939 

-0.2 

8.637 

-3.6 

2.22 

41.25 

9.685 

9.513 

-1.8 

9.010 

-7.0 

1.16 

34.45 

10.17 

10.01 

-1.6 

9.209 

-9.4 

-1.2*  -6.7 


FIMASW3 

4.35 

55.08 

9.379 

9.384 

0.1 

9.223 

-1.7 

2.22 

45.50 

1015 

10.11 

-0.4 

9.813 

-3.3 

1.16 

36.67 

10.73 

11.18 

12 
+  1.3* 

10.62 

-1.0 
-2.0 

FIMASW4 

4.35 

52.83 

9.531 

9.558 

0.3 

9.528 

0.0 

2.22 

44.76 

10.15 

10.22 

0.7 

10.17 

0.2 

1.16 

36.16 

10.67 

11.18 

18 
+  1.9 

11.08 

12 
+  1.5* 

FIMASW5 

4.35 

54.28 

9.405 

9.453 

0.5 

9.367 

0.4 

2.22 

45.39 

10.08 

10.14 

0.6 

9.985 

0.9 

1.16 

36.69 

10.61 

11.20 

16 

+2.2 

10.90 

21 
+  1.3* 
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File  Name 

Flow 
Vel. 

Film 
AT 

Inst 
Tube 
Data 

Fixed 
Method 

% 
error 

Float 
Method 

% 
error 

vw 

(m/s) 

k 

h°2 

(kW/m2-K) 

h0 

h0 

FIMASN4 

4.35 

46.24 

9.947 

10.07 

1.2 

9.512 

-4.4 

2.22 

34.27 

11.12 

11.34 

2.0 

10.24 

-7.9 

1.16 

24.82 

12.21 

._ 

__ 

__ 

__ 

+  1.6*  -6.2 


FIMASN6 

4.35 

45.93 

9.872 

9.823 

0.5 

9.302 

5.8 

2.22 

34.40 

10.84 

10.86 

0.2 

9.898 

8.7 

1.16 

25.52 

11.81 

14.35 

21.5 

+7.4 

11.77 

03 
+4.9* 
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APPENDIX  E.    UNCERTAINTY  ANALYSIS 

In  the  measurement  of  a  physical  quantity,  there  will  always  be  a  difference 
between  the  measured  value  of  the  quantity  and  the  actual  value.  The  magnitude 
of  this  difference  depends  on  the  accuracy  of  the  measuring  device  calibration, 
operator  experience,  environmental  effects,  etc.  Eventhough  the  error  associated 
with  a  single  measurement  may  be  rather  small,  the  error  may  grow  to  substantial 
proportions  when  combined  with  other  measured  quantities  in  a  given  calculation 
scheme.  The  best  estimate  of  the  difference  between  a  calculated  or  measured 
quantity  and  the  actual  value  of  the  quantity  is  known  as  the  uncertainty. 

The  uncertainty  may  be  estimated  by  the  method  of  Kline  and  McClintock 
[Ref.  24].  This  states  that  for  a  quantity  R,  which  is  a  function  of  several  measured 
quantities  (R  =  R(x1,x2,x3,...,x„)),  the  uncertainty  in  R  is  given  by  the  following 
relation: 


•v 


ydx,     ' 


\2    / 


™w2 

ax,  2 


\2 


ax,   3 


/  j 


1/2 


(E.l) 


where: 

WR  =  the  uncertainty  of  the  desired  dependent  variable 
x1,x2,x3,...,x„  =  the  measured  independent  variables 
W1,W2,W3,...,Wn  =  the  uncertainties  in  the  measured  variables. 

Georgiadis  [Ref.  5]  gives  a  complete  description  of  the  uncertainty  analysis  for 
this  experiment.  The  uncertainty  analysis  program  written  by  Mitrou  [Ref.  25]  was 
used  to  calculate  uncertainties.  The  uncertainty  for  runs  using  an  insert  ranged  from 
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±2.2%  to  ±4.6%.  The  uncertainty  for  runs  without  an  insert  ranged  from  ±  10.9% 
to  ±17.1%.  Sample  outputs  of  the  uncertainty  evaluations  are  included  in  this 
appendix. 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS 


File  Name: 
Pressure  Condition: 
Yapcr  Temperature 
Water  Flow  Rate  (%) 
Water  Velocity 
Heat  Flux 

Tube-metal  thermal  conduc 
Sieder-Tate  constant 
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DATA  FOR  THE  UNCERTAINTY  ANALYSIS 


File  Nans:  FNMAUSHB 

Pressure  Condition:    Uacuum 
I'apor  Temperature  ■ 

Water  Flow  Rate  (%) 
Water  Velocity 
Heat  Flux 

Tube-metal  thermal  conduc. 
Sieder-Tate  constant 


76.541  (Deg  C) 

-  20.00 
1  .  1  S  (  m  /  s  ) 
3.9B5E+05  (W/mA2) 

-  385.0  (W/m.K  ) 
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naya  cnp  THE  UNCERTAINTY  ANALYSIS: 
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APPENDIX  F.   INSTRUMENTED  TUBE  CONSTRUCTION 

The  medium  diameter  instrumented  tube  of  Poole  [Ref.  6],  with  six  wall 
thermocouples,  was  fabricated  by  welding  together  a  copper  tube  in  three  pieces 
after  imbedding  capillary  tubes  at  mid  depth  in  the  tube  wall.  Teflon  thermocouples 
were  inserted  into  the  capillary  tubes  for  tube  wall  temperature  measurement. 

An  attempt  was  made  during  this  study  to  construct  instrumented  tubes  of 
small,  medium,  and  large  diameter  in  the  following  steps: 

1.  Take  thick-walled  copper  tube  bar-stock  and  machine  to  specified  inside 
diameter,  outside  diameter,  and  length. 

2.  Cur  four  evenly  spaced  slots,  of  given  depth  and  width  to  accommodate  metal 
sheathed  thermocouples,  from  halfway  along  the  tube  to  the  end. 

3.  Solder  metal  sheathed  thermocouples  into  groves. 

4.  Cur  copper  strips  from  another  tube  to  fit  into  the  top  of  the  slot. 

5.  Clamp  copper  strips  to  the  slots  (using  jubilee  clips)  and  solder  in  place. 

6.  Turn  off  excess  copper  from  strips  to  original  outside  diameter. 

7.  Send  tubes  to  plating  shop  and  plate  the  whole  tube  to  a  given  thickness  of 
plate. 

8.  Return  tube  to  NPS  machine  shop  and  machine  back  down  to  desired  outside 
diameter. 

The  process  was  completed  through  step  2  of  the  above  procedure  and  the 
tubes  are  ready  for  step  3.  Several  of  the  latter  steps  were  attempted  using  a 
practice  tube  with  monel  wire  placed  in  the  slots  on  the  tube  vice  thermocouples;  this 
was  to  ensure  the  process  was  safe,  and  to  prevent  destruction  of  assets  (small, 
medium,  and  large  diameter  tubes  machined  to  specifications,  and  the  associated 
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metal  sheathed  thermocouples).  The  following  summary  documents  the  steps  taken 
and  lessons  learned. 


1.  Completed  tube  fabrication  through  step  2  of  the  fabrication  procedure 
(above).  A  practice  tube  was  used  from  this  point  to  continue  the  procedure. 

2.  Cutting  the  copper  strips  proved  difficult,  so  the  monel  (similar  to  the 
stainless  steel  sheath  on  the  thermocouples)  wires  were  silver  soldered  after 
being  pinned  in  place.  No  lead  was  permitted  in  the  solder  since  copper  will 
not  plate  to  lead,  but  will  plate  to  silver. 

3.  The  practice  tube  was  sent  to  a  local  contractor  for  electroplating.  (Bay 
Custom  Chrome,  Marina,  Ca.) 

4.  Two  different  plating  procedures  were  used  on  the  tube: 

a.  First,  a  cyanide  bath  treatment  was  used  to  electroplate  the  tube. 
However,  this  procedure  severely  scorched  the  surface  of  the  tube, 
resulting  in  the  return  of  the  tube  to  the  NPS  machine  shop  for  repair. 

b.  After  repair,  the  tube  was  again  sent  off  and  subsequently  treated  with 
an  acid  bath  procedure,  which  resulted  in  an  acceptable  tube  surface 
finish. 

5.  The  original  thickness  of  electroplate  was  not  thick  enough  for  our 
application,  so  another  acid  bath  treatment  was  performed  to  bring  the 
surface  thickness  to  the  desired  level. 

6.  The  final  product  was  suitable  for  our  intended  application  with  the  exception 
of  the  following  items  which  need  to  be  corrected: 

a.  Existing  voids  (from  non-uniform  soldering)  or  flaws  in  the  tube  surface 
would  not  fill  in  to  give  a  uniform  thickness  around  the  tube. 

b.  The  practice  wires,  tube  end  pieces,  and  tube  interior  were  not 
protected  properly  during  the  electroplating  process,  resulting  in 
unacceptable  copper  deposition  on  critical  areas. 

It  is  recommended  that  the  instrumented  tube  efforts  be  continued  with 
particular  attention  to  the  following  items: 

1.  Silver  soldering  is  an  acceptable  method  of  fixing  the  thermocouples  in  place, 
since  the  thermocouple  melting  temperature  is  ~-1700°C  (silver  solder  is 
applied  at  ~1100°C). 
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2.     If  the  voids  left  in  previous  soldering  efforts  persist,  then  cut  copper  strips  to 


If  the  voids  left  in  previous  soldering  efforts  persist,  t 
place  over  the  thermocouples  to  hold  them  m  place. 


3.     Ensure  that  critical  areas  on  the  tube   are  well   protected  during  the 
electroplating  process. 
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APPENDIX  G.    DRPINST  PROGRAM  LISTING 

The  program  DRPINST,  which  was  used  to  collect  all  instrumented  data,  is 
listed  in  this  appendix. 
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1(71/171  NJ  =  MnM_  J  MCTpi  IMCMTCn 

1  (7M  Q  I  7D(1                                          TIIDC      TVpp      M  =  MCnTIIM 

1i7l/!CU  ATU                                     DDCCCIipp     rnwnTTiriM 

I    %J  T  vV    .  Til.                                                                         rr.Lk/WWl.L          WWI«l_/XlXWI'« 

10C0I  A  =  ATMnpi-|5p  JQ     (101      kPA) 

10S1 '  U=UACUUM 

1/7.C9I  I.IUCM     A      R       tl      ACpCAR      THCCTUCD      AC       TM 

1   %j  _j  x.    :  m  il.1.       11       u        v        iiiik_nr\        i_/wl_iiil_i          Ilk.        in 

1 053 1  FIMAUSH1    IT    REFERS    TO   A    SERIES 

1054'  OF    RUNS    AT    UARIOUS   PRESSURES    AND 

1  0cn  1  W|&pqp    iJELQCITIES. 

10SE1  5TH                          TUBE    TYPE    S-SMOOTH 
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1058! 
I  059  ! 
!  0E0 ! 
105!  ! 
1052! 
1053' 
1054  ! 
1055! 
1065! 

1IJC7I 

I   xJxj    i 

1  (ICO  I 

(Www; 

1059' 

1  070 ' 

1  07  1  ' 
1  (7i  n->  i 

I    XJ    I    *_     1 

1  073  ' 

1  01*7/1   I 

i   xj   i   -r 

1  (7I-7C  I 

i      1-  ■_ 

1  075  ' 

1  (71-7*7  I 

1078' 
1  079  ' 

1  (7iQ  1    I 

I   \J  XJ   t 

1  {7IQ7 

1-        w       t_ 

1  093 

i  094 

1  0QC 
1  (7IQC 

i    *•  w  ■_• 

1  097 
1  088 
i  0gg 

1  090 

1  (7lQ  1 

t     XJ    *_l     I 

1  092 
i  093 

i  094 
1  095 

1  (TiQC 

(7ia*7 

(TlQQ 

099 
100 

1  (7i  1 

I     \J     I 

1  (717 

I     1/  J- 

1   1  (7i~ 

I     I    -w  w 

1   1  (7i/l 

i    t  \j  -+ 

1    1  f7lC 

I        I     t/  U 

1   1  tftC 

I        1     \J  U 


F-FINNED 

TYPE  OF  INSERT  INSTALLED 

N-NONE 

W-WIRE  WRAP 
(j.u'TATrry 

THE  HEATEX  INSERT  WAS  USED  EXCLUSIVELY 
FOR  FINNED  TUBE  DATA,  SO  POSITIONS 
5  & 7  REFER  TO  THE  TUBE  FIN  SPACING 

05-0. 5mm 

10-1 .0mm 

15-1 .5mm 

20-2. 0mm 

RUN  NUMBER 


CTMAI        DnCTTTnKI 


MCAMTMC      ntr      Al    I        CI    fiCC      TM      DDnCDAM 

I    lb.1    III^IJU         W  I  lll_L.         i      Li    iuj  ±    IX  I       I  \  W  W  t\  M  I    » 


FLUID  TYPE 

OPTION  WITHIN  PROGRAM 

INPUT  MODE 

CTMwpn  nc  CMnnTu 

TMcrrpy    JvpE 

LuUl  I  1  W     «  wXltl-LI.  I        I   \  W   W  I   \  I    it     I 

TIIDC      niAMCTCD 

I    U  i_/ 1_        Uliniu    I    Lm 

ppcccijptr    CONDITION 


TCT 

J.   ■   I 

ISO 
IM : 
IFG 

INN 

TUT 

j.   yv  i 

ITDS 
IPC: 
jnu  ; 
rriM     /rr  /    n  n  \ 

Uui  i       t    xjx*  i         xj  v     i     / 

COM  /Cc55/  T55(5> 

COM  /CcEE ''  T5C(5) 

COM  /Cc57/  T57( 5  ) 

COM  /  C  r  5  9  /  T  5  8  (  5  ) 

COM  /Fid/  Ift 

DIM  Emf(13)  Tw'B)  Tuio(B)  Tw4(c) 

nATA     (7i      1  (7l(7lQCf7lQ  1       7C7?7      Qi^CQ      _*7C*77:/1C      Q7QC      *7Q(7l7CCqC      Ot 

qata  _g7474gccgg  c_g7cggc4.11  -2.6E192E+13  3.94078E+14 
dcah  r t  m.  \ 

1C      7      CQ/17C-7      _*7      7C71C-7     **      7Q/1C-1   1       _Q      771DC_1C      Q      *7171C-7(7( 


*  \ 


nATA    773 

READ    T55 
C^ta   773 

DCAn      TCC 

Mt-IIW  I     UU 

PlATA     777! 
lJ  i  i  i  n    4-  i  w 

READ  T571! 
qata  273 

READ  T58( 

n«=  -p-  i  c7/i    i 

A^»DTtn..P"7 /* 

L- . 1 3335      ! 
LI -.060325    ! 

I  7-   (TiT-zlQ-JC 

PRINTER  IS  1 

DCCD 


1C      7      CC7QC-7      _C      QOC7C-7      -~7.       1  7 /17  C  _  1   1       1       777CC-H      -1       (7)1QDC_1D 

*     ) 
1    C 


Cia7AC_7      _7      7Q77C.7      7      QC7CC-11       1       QTpC-IC     _7      74.QCC_iq 


IC      7      CQ^1C_7      _7      C7^7C_*7      A      (7lCp*7C_11       _1       7  *7  Q  1  p  _  1  C     C      /  /I  Ol7  C_  7  (71 
*    ) 

'  Inside  diameter  of  stainless  steel  test  section 


!  Condensing  length 

!  Inlet  end  "fin  length" 

'  Outlet  end  "fin  length" 
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i  07 
1  aa 

i    \j  \J 

1  (hQ 

1  (A 

■   w 

1 
1  2 

7. 

4 

"7 

20 

23 
25 
29 
1  32 
1  35 
36 

A   1 

44 
47 
52 
53 
cc 

7  1 
74 
77 

C7 

99 

QC 

98 
20  1 

4-  \s  -r 

207 

7  1  I7i 


2  I  I 

2  1  3 
?  1  p 

2  1  9 

220 
o?  1 

o?c 

225 
227 
228 

990 


t  mdi  it    "ccirrT    ci  inn    /ra.cTcaM       i=p_i  17       ^.pD  \ "     t  *  + 

DDTMT  IICTMC   "  A  V   "C.l.^t   nn+inn*""" 

1  Take  data  or  re-process  previou: 


PRINT  USING  "SX,""  i  iars  coxa  or  re-process  previous  aata 
PRINT  USING  "EX,M"  2  Print  rsu    data  (NOT  COMPLETE) 

TMDI  IT      T  «~ 

TD      T-«\1       TUCM      7fflm 

A  I  A   O  W  .-      I  lMi_ii        w  W  w  i_ 

PRINTER  IS  701 
Ijob-0 

PI  DAP   "7(7iQ 

BEEP 

tmdijt    " dmtdd    fvinMTu      DftTE    AND    TIME    ( MM : DD : HH : MM : SS  ) "    DateS 

niiTDiiT    "7(7>a  .  "  jn  ■•  .  n»  +  «* 

ni  itdi  it    -7171a  .  "  Tn  " 

W   W      I     I        W     I  I     U    J     1  I     L-f 

TD      T , ~K= 1       TUCM 

A   ,  A    J    W  i->  I  llll—lJ 

BEEP 

TMDIIT      "CU'TD      DACD      AMP     UTT      pMTDD"      f\L 
DMn      TD 


dmtdd  709-DsteS 

DDTMT   " 

PRINT 

DDTMT  II C  T  N|  0   "10V 
TD   T 

DDCD 


MnTD 


«_icC   ^-iiii^t    i..LrK 


nar 


.  n  «  +  -  * 
DRPINST 


TMDIJT   "DMTDD   TMDIJT  MOnC   (  I  «M[U  OAfA    2=EXISTINS  PILE)"   I  ? 

BEEP 

TD   T»,«1   TUDM 

BEEP 

TMDIJT      "gTUD      A      MAMC      DQD      TUD      DAM     QftTA     p  T  l_  D  "      n      f^»* 

PRINT  USING  " 1  EX  ""Pile  nane     :  ""  14fV';D  file* 

QDDATD   gnAT  n   f , 1 £ J   ID 

ASSIGN  @Piie  TO  0  file$ 
BEEP 

TMDIJT   "DMTDD   DDQMDTpy   CODE   '  !  "F I NNED  0~PLftIN)"   Ifc 
pDTMTDp   TC   1 
T  „  ^.  =  f7i 

BEEP 
PR  INT 
PRINT 
PR  INT 

DDTMT      "         't 
DDTMT      " 

TMDIIT      T-.-, 

a  i  «i    w  t       a  i  1 1  i 

OUTPUT      AC, l..Tf„      Tn„ 

tmdijt  "mq_  np  THERMOCOUPLES  IN  WALL"  Inu;t 

TD   Tfr.l   TUDM 

INPUT  "ENTER  FIN  HEIGHT,  FIN  PITCH  AND  FIN  WIDTH  (Fh,Fp,Fw>  IN  MM",Fh,Fp,F 


SELECT  INSERT  TYPE:" 

(7i  MflMC      I  nCDAI  II    T  \  " 

w        iJ Jnu      v  ,_/ 1_ ,   n ^ i_  ,    / 


1 


TWISTED  TAPE 

WTDD  WD  AD" 

HEftTEX " 


0  30 

9-  1 


ELSE 
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'  «.<3«-  I  H  ^ 

1  9^/t  c,..ra 

1235  END  IF 

1237  OUTPUT  ©FileiInwt,Fp,Fw,Fh 

1 243  ELSE 

1245  IF  Ijob-1  THEN  1255 

\  r>  aq  ocpp 

i  4.  -T-  *j  wutir 

1252  INPUT  "GIUE  THE  NAME  OF  THE  EXISTING  DATA  FILE"  ,D_f i le$ 

1255  PRINT  USING  "1EX,""This  analysis  was  performed  for  data  in  file  "",10A";D_ 

f  lieS 

1259  IF  Ijob-1  THEN  1270 

251  Nrun-18 

2£4  BEEP 

257  INPUT  "ENTER  NUMBER  OF  DATA  SETS  ( DEF-1 8 )"  ,Nrun 

770  accjckj  @Fiie  TO  D  file$ 

282  END  IF 

7>OC  TC      Tinh-1      TUCM      IU7 

9QC  DDTMTCD      TC       1 

l-  ~/U  r  ~\  1  M  I  Li  i       a  J        i 

340  BEEP 

3c ]  PRINT  "SELECT  TUBE  DIA  TYPE:" 

7C7  DDTMT      "1  CMAI    I     " 

«/^ji  1    1  \  x  1  ^  1  t  ^1  i'  ill 

3gg  ppjKiT    "2      MEDIUM    (DEFAULT)" 

aCa  pOTKJT      "3  J_ApgC" 

ICC  TMDIIT      T4-H- 

3.CQ  pDTMTCp      TC      701 

375  Di=.0127      '  ID  of  medium  and  lar^-e  tub°s 

■no  n ~ «    ffliapc  1     nn    ~f    ».  •  h  <  i  •  •»    +  . .  ^  ~ 

yi  c  a  ni=    iTiiarac  t     nn    ^<=    ,  .^l,,.,,  +  »^i     ,„i^  +     i«^r,+  ^ 

T  — '    -J  L^     1  ■    U    I      JVJ  W  L-1  w    i  Ul   M   I^W     b    WU  J.    ■    *  .t    W    I.  X    w  I    I  Lj    *<    I    I 

/ict>  n?.    (income  1     nn    ~£    , ,nhi»^+i»H    r-"+'i»+     i-^n+h 

—  »_'—'■  i_  .■__■..     i     _.  '    ■  —  ■      riWwwL4|_rj.w        uUi      xui         UCp    d  i  i 

455  IF    ltd"")    THEN 
acq  n  ,  .    aaacitz 

~T  w  lj  ui  «    t_  w    _    _■  *-    — 

TU     I  UU  ■     %J     l     i_      1 

472     Or-. 01       !  TO  BE  MODIFIED  WHEN  KNOWN 

474  p  Ki  n  t  p 

1  477  jp  itds-3  THEN 
1 47S     Dc~ . 025 
1  479     Dr*"  .  02 

1400  END  IF 

1  A  QC  l/^,.«70C 


'TO   BE    MODIFIED   WHEN   KNOWN 


1498  Rr»!-Do*LOG(Do/Di  >/<2*Kcu>    !    Wall    resistance   based   on    outside    area 

1  Cfll  1  DCCD 

'     -1  \J     1  CJ  L-  L_  1 

1504  INPUT  "ENTER  PRESSURE  CONDITION  < 0-U  ,  1 -A  ) "  ,  Ipc 

1510  BEEP 

1543  PRINT  USING  "1GX,""This  analysis  includes  end-fin  effect 

1545  PRINT  USING  " 1 5X  , " " Thernai  conductivity      -  "",3D.D,""  (W/rc.K) ;Kcu 

1549  PRINT  USING  " 1 SX  , " " Ins i de  dianeter,  Di       -  "".DO.DD,""  (mm) ;Di*1000 


111 


1  cc? 

1  can 

icua 
■  *_/  u  *j 

1  C7C 

1  £39 
i  g42 

1  c  a  a 

t     t_/   -T  \-/ 

1  CC  1 

i  C54 
i  557 

1  CCA 

w  w  «. 

1  CC7 

1CCO 

i  572 

i  575 
i  575 

1  CO  1 

2  .  DD ' 
534 

7  0*7 

77 1 
774 
777 

-7001 

i    w  W 

763 

734 

735 

797 
70  a 

792 

795 

799 
ara  1 

W  W     I 

9r7l4 

W  U/      I 

3  1  2 

9 1  3 
c  1  c 

W     I     w 

9 1  9 
922 

CJC 

W   1-    wJ 

- 1  - 

93  1 

Q7J 

W  ^  T 

640 
943 
849 
q/i  a 

w  -r  u 

9C2 

QCC 


PRINT  USING  "16X, ""Outside  dianeter,  Do      -  "".DD.DD 

IF  Ijob-1  THEN  1648 

BEEP 

INPUT  "GIUE  A  NAME  FOR  WALL  TEMPERATURE  FILE",Wtf$ 

CREATE  BOAT  Wtf$,E 

ACCTCM     «C.1«1       Tfl      U  +  +~<t 

.■-_■■  \Z  I        J-      -i     —      '  I     W  "v     1.      1      — 

TC      T   ,  »U_  1       TUCM 

At  iJWU  I  I    iil.ll 

T   ~  w  =  7 

GOTO    1702 

END    IF 

BEEP 

tmdijt    "SELECT    OUTPUT    (  1  =SHORT   2=!_0N6    ^=RAU   DATA)"    lev 

TC      T  <>„«0I      TUCM 

PPjmt    USING    "  1  5X 

CI    cc 


^n 


;Do* 1 000 


t,  ,^«»    +  wr« 


■       TMCTDl  IMCMTCn      CMnnTU"" 


T,  ,U  •»      +  .,,-,« 


TMCTpjjMC^Tcn    FINNED"" 


CD  TMT      I  IC  TMC      "   t  C  V 

1     I  \  A  1  *    i         w  ^  *  i  *  w  iuf' 

DDTMT     IICTMC      "ICY      "  "  C  -  ^     ~  ,  +  r-  1-,         i.iiH+h         =  n A     U->.^^  +      I  ^^  \ 

.Cn        C.  .       CU 

1 1   H  ,  1  >~  ,  I   ' ' 
CNin    tc 

J-0 

D«»r,~  =  +     ■      I 
'»  w  )J  C  u    1.    • 

T  .  J  4.  1 

ye    lr,=  ]    THEN 
BEEP 

nilTDIIT      "7f7iQ."AD      (\C  AO\      Al    A  (71     tlDC" 

WWII       W     I  IWWJ  I     I  t   \  III        'TU  I     1  l_  T  W  V    I   \U 

OUTPUT    709; "A9    9A" 

TMDIJT  "  l_  T  L.  ;C  TQ  Cf-jCQL;  we  CONCENTRATION  (  1  =  Y  0"=N)?"  N" 

IF  J = 1  THEN 

nijTpijT    70g."Ap    AC42    ft[_41     UR5" 

END    IF 

DCCD 


nn    nn 


nn    iv    7    nn    ?v 


TMDIJT   "CMTCD  FLQUMCTCp  DC  An  TMC"   C  r 

OUTPUT    709" "AR    AFE0   ALE2    UR5" 

ni  ITDI  IT      7f7iQ."AC      C^" 
CMTCD      "7(710  .  C  +  r, 
QMTDIIT      70Q."AC      CA" 
DCCD 

tmdijt  "CONNECT  UOLTAGE  LINE"  Or 
c  m  t  c  p  7  0  9  ■  B  v  o  1 

BEEP 

TMDIIT      "nTCrnMMCPT      lini    TACC      I     TMC"      nt 

Alii      w    1  l-a^wwiiiiCwI  x/wi—iMw*—        u.  a  i  J  l_         |^'s 

TC      P„nl  ,        1       TUCM 

J.    I  W     V     W    A       V      •        I  I      I     I  t_   I    » 

BEEP 

tmdijt    "TMiiAi  tq    U0LTAGE    -    TRY    AGAIN'  "    OF 

GOTO    1819 

END    IF 

I       Q'JTPUT      7gQ."AC      CA" 
I       CMTCD      TGld  .  11+ -»r. 
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1  DCQ 

t   w  u  w 

1  OC  1 

1962 
1867 

W    I     «J 

D-7C 

W     I     Uf 

879 
882 
9°5 

QQQ 

^  ^  *_' 
89 1 

QQ/t 
w  <-i  ~r 

897 

Q(J7. 

-'     \-.       w 

905 
91  2 

91  5 

a  i  c 
w  i  t-i 

91  8 

Q-7  1 

*_»  i_      I 

924 

Q"!7 

930 
g-7- 

935 

939 

942 

945 

949 

95  1 

954 
acm 

ac  1 

JU     I 

952 
972 
984 

QOT 
3  w  t 

990 
99* 
995 
999 

(TifTi? 

*>  w  c- 

-PfTifflQ 

_    I-1    <-  ■     _< 

2009 

201  1 

ora  1  -> 

—  »- >  i- 

r>oi\  a  i 

l_  V   I   *T    ! 

ora  1  c 


OUTPUT   709; "ftS   SA" 

ENTER  709;Bonp 

Etp-Etp*1 .E+6 

OUTPUT  709; "AR  ftF40  AL53  UR5 " 

Nn-7+Inwt 

FOR  1-0  TO  Nn 

OUTPUT  709; "ftS  Sft" 

IF  I>7  THEN 

Sc-0 

FOR  K-1  TO  20 

ENTER  709 ;E 

Se-Se+E 

NEXT  K 

ELSE 

ckitcd    •yaa.c 

C»,f  (   T   ).ADC(  C  ) 

END      IF 

C«fMWC«f/|   \*l      CiC 

l_lt   I     \    X     /        Li'i   I     v    1     /    •'     I     *   i_    '    D 

NEXT    I 

nilTDIIT  7UQ.  "AC  CA" 
m  ITDI  IT  7I7."TI  D9C  " 
I.IATT     r> 

W  t  l  A   I 

ENTEf 

OUTPUT  7  1  3  •  "  T^R71- " 

MATT   T> 

Wl  I  X     I 

ENTEi 

OUTPUT  713* "T 

UATT  ? 

CMTCD  7]~;Ti2 
Ti=(Ti  lijn  \ * 
p 

BEEP 

TWDIIT      "CMTCB      DDCCCIIDC      CAI  ICC      DCAnTMC      /Dn,   \"      D^» 

A  I'll       W     I  l_  l  1     i     _  I  ^      _     _j    w  i       ^  W  it  *J  W  l_  lUllJlllU  II       y  G      /  -i    - 

D„.nl«Dn.»COQ/l      "7  IDCT      TO     D  = 

t       ■-     -  ■  f.  ■     ■  w  w  ■„    -     _  ■  -r     .  :    I      w/  A  I    w        I      □ 

ELSE 

CI\ITCp     «c^ 

tc     t« 1     no 


,  an,i_i.hi,i_rii  v      /,■  n.i  ,  ,.i_,i   >uH 


r. 


9£7l      nD       I«M^,.r,      TUCM 


ELSE 
"a     ^ 

END  IF 
END  IF 
Ts  t  csr;  1  "FNTvsv57(  E~f  (  0  )  ) 

t  «  +  - a ~ 2  =  F N T  v  s  v 55  (  E m  f  (  1  )  ) 
Tst2ar,2-Tsteam2-273.15 
Tateam-(Tateaml+Tatean2  )/2. 
Tstcan-Tsteanl 

T^nni».rKIT„«wCQ/  CT-fC  "7   \   \ 

IIWWI'I  I         I    *      1        t       _'      V      ^J   *_■      I       ,_    I     ,      I  .       4_         •         - 


113 


20 1  S 

20 1  9 

t     -.  <—  ^ 

202  1 
20°2 
2023 
2024 
2026 

20™2 
7044 
2057 
2059 
2059 
2050 
2052 
-?055 

_.  u   — i  _■ 

k-UU      ■ 

?(J7  1 

2072  ' 

20?3 

2075 

X.   fc/     ■        I 

2090  ' 

205^ 

-  0!  2  9 

1 0  9  3 
7/7104 

209° 
!  0 

1     A 

1  c 


Tccrs~FNTvs v5S( Emf < 7 


!  r  c  o  n  -  z  ,' 


-nn-z    1  c 


1  c 


40 

/A 

A  C 

-r  «j 

/t  Q 

50 
52 


Tccn=Tcor!-27L. 

T^^-g  _ 

crno    t  „(7i    m    Tm.i+  -1 

T...I    T    \.CMT„.,X7/  Pff  '   T-»-0   \   \ 

toj(  t  )»TW(  t  )-273.  15 
Twrc-Tww+Tw*  I  ) 

MDVT   T 

I  «1_  /V   I       X 

D  =  .  +  -PMC  ..  =  +  ,  T«  +  ..~  . 

OUTPUT  709" "AR  AFE4  AL54  UR5" 

nilTDIIT      "7(710.   "AC      CA"  I       DQPCCIIDP      TDAMCniipCD 

Ss*0 

pnp    i/  =  i     tq    20 

ctiijrp    709'Eirsri 

Ss»Ss+Etran 

NEXT  K 

p+r  =  r,  =  Agc(  s=/20  ) 

BEEP 

ppcccnpr    im    p3    ppnjvj    TRANSDUCER 

pv3^2  =  ( -2 . 93S04*Ptrarc+ 1 4 . 7827  S*S994 . 7 

i  Gft!jsp 

1  TRANSDUCER 
j  cnp  ft IR  CONTAMINATION 
Pkn*Pvftn?*1  F-3   I  TRANSDUCER  IN  kPa 

I       CAT  DDPCCIIDP      TM      LDa 


is  P      .     /up 


FF 

T  «= 


T,,=r-,"?=:PMTw=rW     D  w   , 

1     »  up  1.       I     n   1    v   J(J  v    1       »  ( 


CAtlCP      TM      L  .p, 


1 ' «  +  *"  c  N l '  v  * 

M.M.elO       (71 


t  <  T  5 1  e 

c 


OIV.       ) 


TP      T  f  +  =  1       TUCM      M,.,„»IQ7      aD 

XI  XI        W  I  I     I     I  [_  i   V  I     IMT    V  r     w     1       .     _'  w 

TP      Tf*=?      TUPW      M,,>.«0      (7IP 

XI  X      I       <-  *-  1     I     I  1_   I   N  I    IUI    V  U   t.     •     t/U 

*J  CJ    i.      '   •      \    1       v   <_*  w     1  \      1      •  \    1    I*. 


/?D      QC   \   \*D, 


TIT      T~w«"3      TUPM 

xi  x  w  v         x.  1    i   1 1_  1  0 

PRINT 
PRINT 

DP  TMT      I  IC  TKIC      "  1  (7lV       "  ' 

i      i\Xl«i  \_/  w   x   i  i  ^/  i    i/A     . 


,4-      W   .'    D  w  =  .-,?-  /      1  _   f    M ,  . ,  , 

U     /  /      U       V   C3  |-i  X.  vi.  iw 


:S.9E  )  >*Psat  ) 


n. 


,et 


jnbe'" 


nilTDIIT      "7(7iQ."  AD      AC  A  ft      Al    A  (7|      HDD 
nilTDIIT      "7f7lQ."AC      CA" 

WWIIWI       I     t^  *-»  1     I  I  »-/     -^'i 

END  IF 
PR  INT 

TP   T~..«?   AMI-!   Kir-..  1   TUPM 

■  i  i  C  i  i 

Psst 


DDTMT   1ICTMC 


1  Y 


DDTMT   liCTMC 


PRIh. 


'3  ) 


1  I  D» 

\  h.  t  wt 


Ptrsr, 


(kPa  ) 


T+  K-  =  ^ 


sns 


'  r  ) 


(C) 


NG  1  % 


Mclal 


NG2%" 
Molal 
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2161   PRINT  USING  " 1 X ,3( 3D . DD  ,2X  )  ,2( 3D . DD ,2X  )  ,2 ( M3D . D  ,2X  ) " ;Pka  ,Pkg  ,Pkp  ,Tvap2 

i       Mfpp  1        Mt^^~> 

t  ,m  ng  ,  ,nn  iyc 

7  1C/1  DDT  KIT 


,Tsa 


2167! 

"3  1  7(7. 


'Energize  the  vacuum  system  


END  IF 

IF  Mfr,g1>.5  THEN 

2173  BEEP 

217E  IF  In-1  AND  Ng-1  THEN 

7  1  -7Q  DDDD 

2122  PRINT 

7  1  QC  DDTMT  IICTMC   "1fl»V 

4_     I    h.   —  I     I\1I1    I  -    J   J.   '   -   ■_  («U/\ 

9  1  DO  DDDD 

t-     I  W  W  L^C<—  I 

2  l  91  INPUT  "OK  TO  ACCEPT  THIS  RUN  ( 1 »Y  0=N>?"  Ok 

2194  IF  Ok-0  THEN 

2 1 97  BEEP 

T?an\  ptcd  -motd-  tutc  naTA  cct  i.iti  i  pc  nTcr-Apnpni 

2203  WAIT  5 

77<7|C  CnTfl      1  -70(71 

—  —  v-  U  U     I   I  \_f  <W 

2209  END  IF 

2212  END  IF 

22 1 5  END  IF 

2219  IF  In" 1  THEN 

7771  TD     Cv./    1  pi     flD      p~\1fllfll      TUCM 

777/1  T  f  ~>«0I 

2227  BEEP 

2230  INPUT  "INCORRECT  FN  ( ! "ACCEPT 

2233  IF  Ifrc""1  THEN  1304 

223E  END  IF 

2239  END  IF 

2242 '  ANALYSIS  BEGINS 

2245  T  i  1  =  F  N  T  v  =  v  5  9  (  E  n  f  (  3  )  ) 

T 1 2"FNTvs v55( Enf ( 5  )  ) 

jc  1  -pjJTvsv5Q(  £~f  (  4  )  ) 


,nci  dtd/  nccam  t  \  \ "     t  £~ 

WU.I— Dtl—   VL/DJ       ilUU.     I       /      /  ,    1     I     ™l 


7245 
2247 

77 /Id 

i_  w  -r  w 

2249 
2250 
225 1 
2252 
2254 
2255 

77CC 

W   A_   _'   — 

77D7 
7  7DQ 


Tn7»CmT„, 


xc(C«{/cn 
—  u   ^w/ 


T  i    1 


1  _777*   1  C 


Ti2"T:2_273 

T~  1  »T~  1  .77"; 

Tc2~Tc2-273 

T^«  1   1  «T~  1  _T 

I    uC  A     r  I    w    i  I 

TrV  1  7=T~7_T 

TH-  1  ^«T?-T  1 

C 


1  c 

15 

1  5 
i 

7 


D+  „  1 


i    v  3  i  ■  k.  i.  P  /   t.  flu  , 


n+W«-7      CQ71D_7_1       DfTl/t  C  /JD_C#D  +  r,  1  +  1       7  1  7(7l  1  D  _  1  r7l  *  D  +  r.  1   A  7  _D      t   1  C  >ip  _  1  C  *  D  +  r,  1   "  1 J.^      7  7(711 


D_1  Q  *  D  +  -,  1   *  A 

77D7  T~T«=T ,   1  xTr  i   . 

77C7         DDTMT      IICTMC      "IV       ' 

*-*_l-JW  llVllHI  UwxihU  I    f\     , 

DELT3   TPI'  E   " "  " 
22E4   PRINT  USING  "IX  ' 


TTM1  TnilTI  TTM7 


T0UT2 


IN3 


T0UT3    DELT1   DELT2 


(T,(U^^ 


(Met  a 

77CC  DDTMT      IICTMC      "7V       1  (h  I  7  n     nn      7V  V'lTil       T^l       Ti' 

—  ~    -   -  I     Mill  I         Jt/itflJ  <_/\}rflU\A.t-/at_rk/|A-/t/  tlli]i\-/i)lxi 

22SE   Er 1 =ABS( T 1 1 -T  1  ) 


l  i 


Jo; 


(Quartz)  

,T1  ,T2  .Tdell  ,Tdel2  ,TdeI3, 


Tns 


115 


2257  Er3~ABS( Ti2-T !  ) 

2258  PRINTER    IS    1 

2272  IF    Erl>.5    OR    Er3>.5   AND    In=1    THEN 

2271  BEEP 

2272  PRINT    "nCT    AND    TC    DIFFER   BY    MORE    THAN    0.5    C" 
2275  Ok  1 " 1 

2278  END    IF 

11Q~)  TC      nt  1  .01      AMn      Cr.1-,      C      AMP      T~«1       TUCM      I^Qf?) 

■79QD  TC      nM.ffl     AMn      C  r1\      C      Min      T«dI       TUCM      17Q0I 

*-     l-      w    1.  Ll        J.  ME  W    \       \        I      tm  11/  I  X     O       /        /     /         \         I      X.  II* 

2293  IF    Er2>.05    AND    In"!    THEN 

2295  BEEP 

??Qa  DDTMT   "HPT   AMH  T-DTI  C   nTCCCD  OV   MODC   TUMI  CV" 

x.  x-  »j  *j  i     i\xiii  u/wi         i  1 1  i  L_  i        i      x  L_  I—       Liii     i     l_i\       ky  I        iiumL         i  null       xj  /o 

2302  Ok  2-1 

-?7f7ic  tc    nii^n    Ai,m    Cr-?\    01c    AMn    t~=i    tucm    1 -7Q0i 

.-     -      v      -  xl  UIVL  "X/         I    II  1  W  l_  I       X_    .■        .     x     -^  I    IMU  Xl'l  I  i     i  i_  i   1  I       IWW 

230S  END    IF 

23-11  D  D  T  jv  I 1 T  C  D      T  C      701 

2312  INPUT    "COOLANT    TEMP.    RISE    MEAS.    ( 0  =  TEFLON      1 "METAL    S.      2=QUARTZ      3=TPILE>" 
1 1  n 

2314  IF    Itn-0    THEN 

2315  ppjmt   USING    "IX    ""USING    SINGLE    TEFLON    THERMOCOUPLE    FOR   COOLANT    TEMP.    RISE" 

2317  T 1 1  =  T 1 1 

2320  T  2  ~  =  T  c  1 

2323  END    IF 

2325  IF    Itr,"1    THEN 

"327  PRINT    USING    "IX    ""USING    SINGLE    METAL    SHEATHED    THERMOCOUPLE    FOR    COOLANT    TEM 

P .    RISE' 

?7?a  T1  i.Tp 

-  ~-.a  CMn    tc 

9330  tc    itro  =  2    THEN 

~>~~G  DDTMT      IICTMC.      "1V       ""IICTMC     DilADT?      TUCDMCMCTCD     COD      0(1^1    AMT     TCMD  DTCC""" 

*_  ^  _/    xJ  l       ivxill  X/  x/   X   I  *  L/  I     I'.      .  \_/WXl1W  Xlf  W  11  i    v     I     X.  llll_l\flwill—     I     UK         I      \_'  I   \         UUuLniH    I  1     L_  J    I  I        »  IVXw/L_ 

2341  T 1 : " T 1 

"'344  "r^o^T"' 

2347  END    IF 

2349  tc    t+m«3    THEN 

7-/IQ  DDTMT      IICTMC      "1V       ""IICTMC     TUCDMODTI    C      COD      mfll    AMT      TpMR  P  T  C  F  ' 

X.   W  "*  W  I      I  \  i.  I  M    I  WwXilW  '    '*     7  Ww'XIIV  tllClMIWI       iuk.        I      Wl\         uu  JuMI^    I  I    Ci   ll        .  hlwL 

2350  T 1  :  "T  ;  1 

T7C  1  T">^  =  T^- 

2352  END    IF 

2353  Tsv""( T 1 i +T2c  >*  .  5 

2354  I f t "0 

?-cc  r„,,«cMrr,,,iT.„n\ 

'■■"'"  v/(j«»i   ■  nwpw\  iuvy  i  < 

77CQ  DU^,..  =  CMDU0.,,/  T»>,n  \ 

._  w   _'   _'  i  tnuuu       i     [mil  Jui  i     I   u  vy   l 

?7C7  M^«,'C      Ti(JQ*(:i.xl7      010  "7  \  /  1  0,0101 

?*CC  M-4  =  I«MW1       0I^CC_1       QCC  /1/(C_^*T  1   ,  J.C      ^CJC-C-TI   .   A?   W  I      PlP^"? 

77CQ  M*-NM  /Dk0,,, 

*_    w   -^     —  I     I    I  I     IU'      I    M    IWW 

*>~zn  1  IK..MP  //  DT*H,  A7  //I    \ 

i_wii  vim     III  /    »  r  a"w*      t>  /  *r  y 

7"*Q^  TC      Tnn.f7l     AMn      U...S      C      TUCM     T')n  =  T?.-/  --)      77C_^4.1       7CC  -A  *  I  1.  .,  +  Q     ^CC-^IL,^-  1       DCC_C* 

*_    w'    ■_■■   w  XI  At   ill  "*_-  '     M    *  t_-  VUU/       •     J  I     lILK  i      —     - '  I     i.  \-  \  i_      .       i      _'   i—  -f       ■        i       .        i      _l  l_  ~T  - ■     **J      '       ^     •     w   w'  i —  *T'VW  i—  i       .      w'  i_l    k_  — ' 
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lJu;"3  ) 

2389  IF    Inn-1    THEN    T2o-T2o-(  -B  .  44E-5+1  .  71  E-3#lJui  +  4  .  45E-4*UiuA2  +  4  .  07E-5*VwA3  > 

2407  IF    Inn-2   THEN   T2o-T2o-( -3. 99E-4+2.75E-3*Vw+1 . 45E-3*0wA2+8. 16E-5*VwA3  ) 

2408  IF    Inn-3    THEN    T2o-T2c-(  8  .  S7E-5+  1  .  23E-3*lJw+  1  .  03E-3*'JwA2  +  8  .  1  6E-S*lv>wA3  ) 
2413  Q-Md*Cpw*(T2o-T1i ) 

2416  Qp-Q/(PI*Do*L) 

2417!  ITERATE  TO  FIND  OUTER  WALL  TEMPERATURES 

2418  Tumi -Tun 

2419  Tuim2-Twn 

2420  Tuim1-(Tuin1+Tun2  )/2. 

2421  Ct-Q*LOG(Do/Dr  >/(2 .*PI*Kcu#L  ) 

2422  Twm2-Twrt+Ct 

2423  IF    I°iBS<Twi*i2-Tuivi1  >>  .001    THEN    GOTO   2420 
2424'  ITERATE    TO    FIND    INNER    WALL    TEMPERATURES 

?   ^    *)C  T  iiMA  K  T  l, IMS 

_    1    —     _  I     1AM    I  W  luJll 

->A->C  TlibsA-T'II" 

a- -t  a.  w  I    iwri-T  <    (Mil 

"?A7I*7  T.   lw,T.=  f    T|.,-7J.T,  ,rY>A    \    l~> 

2428  Cti**n*LOS(Dr/Di  )/(2.*PI*Kcu*') 

2429  TLjn4=,Tu;n-Ct  i 

t.  -r  w  w  j.1  I  lu  u  iuuiiS-  I  UinJ    //     .    v.   w    I  lili—it       W  w,    w        A.  -T  A.    I 

oat.\  tv.»„»0i 

2432  Tu;rii=0. 

2433  FOR  I="0  TO  Tnw4'  — 1 

2434  Tuic(I)",Ti*!(I)+pt 

2435  Twi(I)"Tw(I)-Cti 

A-  -T  W  W  I    Wl'lU  1    OVI'IW    *       I    WWW    *     1      / 

2437  Tu;ni  =  Tujivi:-l-T^i(I) 

2438  NEXT  I 

r  A~a  TW!v.c=T..ify1«/TpW  + 

?  4  4  0  T- ■  •  ■  ■*•  i=Tujri'/Tn,:t 

244]  pDTMT 

2443  PRINT  USING  "IX  ""Position  riuitber  :  1      2      3      4 

5      E  "  "  " 

•7  ,i/i/i  ddtmt    iictmc    "iv     »  "i  «„■  l     „,.+  ..-    ,,.n     +<»^n«        /  ncc    pi        ■     " "    c  /  nn    np     ivi".T,.r(« 

A-  ~r  T  -r  ,      ,  \  a  i  *    ,         w  w  A  i  s  w  ,    n     ,  L_  w  w  w  A         w  w.    ,.    *-  ,  vj  O  A    A  ».   w  ,   i  |J  w    .  '    w  L_  w        w    /  ■  ,U    \    Uu    ■   uL    ,     '    '^     ,  1     I    W  w   v 

\ 

2445  ddtmt 

OAA1  DDTMT      IICTMC      "IV       ",,flw«r.»n-      n.,  +  »r      =.  n  H      ,„„«,.-      ...  -  ]   1       +«~r,«=      _      "  "      T/nfl     np      1  V  V  •  Tfill"l 

A-  -r  -r    .  ,     ■  v  a  i  *    ,         w  w  a  i  .J  w  i  A    j  t.vw,     w  w,  w       w  w    -    w  ,  w  I  i  U        ii  iiiu,  lv  G  a   a         w  w  i  i  aj  w  }A.\wUaU/w,iA/         1    I   wn 

~       T...M, 

W      j      .     WW.    I   J. 

7  4  4  g  Dp  T  KIT 

■?  A  C(7i  T  *  +  «r7l 


a 


245  1 

2453  ML.wi™FNMuw' Twrci  ) 

2454  Rsi«:Rhcu;*ULd*Di/M,'u) 

2457  Prw-FNPrw(Tavg) 

2458  Dtui-Q*LOG(Dc/Di  >/(  2*PI*Kcu*L  >*  .5 

O  ACQ  I    wH«(T7n-T1   ■    Wl    flC/   (T,,»,-Tt   ,    \/IT,,«,-T')„\l 

-—-_'_  i_  n  v  w       V     >   a.  w         iia//l_ww\viUmiI         ,     i*//\1Uuiia  I  i.J    r    / 

?/ICf7l  D.^...DT.n, 

ATUi;  i      Ci       All        I        .  ui 

7/1C1  C,.^  f  .  ,  .DT  *n  i  *l 

A.  T  W     <  WW,        ,     W   A  ,        A^WX^i— 

bTUb    !  J-  '  *    ,         ,\Ca,..Ivjw'aU/?u-,ic-WjiC      ir,    ,  aU,       .    G  A     ,  Uf 
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2457 

2456 

2459 

2470 

247  1 

2472  ! 

2474 

2475 

2475 

2477 

2478 

2479 

2490 

?«ci 

2482 

24S3 

9/10/1 

*_  *T  W  -T 

2495 

2495 
■?  a  an 


2499 
2499 
2490 

i.  JVt' 

252  1 

9C3"7 

2  5  3  0 

~  r  "E 

2575 
2576 

2557 
2590 
2553 
2595 
2599 

TCffl? 

2505 

«_  ^  I.-  w 
"?C  1  1 

*.  O  t   t 

?C  1  A 

2524 
2634 
2535 

...  -j  -  •-> 
->C~Q 


An.r.,1«DI»inr"'7-n,    "7    \/^ 
I   ii   lui       w   ■„      I  t  l.     » ,  *_         l_     i  —     /   *     *T   i 

i  i  i  i  l_x  i     C  w  4-       '  L  '  i-       LJ  x        —  -    - 

u,-n/cc..^f-i*!  «i^i^ 

I        I       I  \    1%  W  Ll  V  J       W  '        j     i     i  i     i   :    i  t_.  i       ! '      ■        '  ■      ■-. 

Ff2-(Kcu*Perirc*Anerea2  )'N  .5 
PRINT    ftnareal  ,Ansres2  ,Hi  ,Ff 1  ,Ff2 
Hil-Hi 

U  .    ■?=(-(  i     } 

ai»i«1.//U,T»p.ri1.\//l,^,,»4r.r..l      \\"         D 

p f  3— p f  1  *Hi  2 "  .  5*  (  Twrci-T 1  i  ) 
p t a=c f2*Hi2"  .  5*  (  Tu;ni-T2c  ) 

C.  .^-C=C  f  7.CMT  =r,k  ^   A  1    =  ~.  1    *1      1      \ 

i   wii^\_/     i    I  w^rin  Oiiiivnioni  Tui    / 

Func7"Ff  4*FNTanh(  ft '  an',*L2  ) 


c  ,^^c  =  u 


*5urfai 


utn 


Funcx  =  FL;nc6- 


inc 


.p.  ,„,-Q_n 


nf,,nr.7Bi 


c. . 


c  /u  • 


)  J.  7  _  »  C  t  a  *  |_  1    /  (    i 

)  +  7  .  *  F  f  4  * '  ' / (  ! 


Cn^^r7  /I 


■FNCo: 


h(  2 


♦ftlanl' 
*ftlar»2- 


LI)) 

i  r>  \  \ 


n  f  l;  n  c 

u ,  1  =u 


DfLir;c5  +  DfL;r;c7+DfL;nc9 
_  (  p  [j  »,  c  .  /  D  f  Li  n  C  X  ) 

tc    adc  /  ui  i  _Lia 2  )>  ,  05    THEN    GOTO    2475 
u ± =lu 1 

DDT  NiT      u  . 
■    i  \ a  1 1  i       'ix 

Ufn  =  CMUfnl     Te   +   «»«*     \ 

i  i  i  y     i   iMi  1 1  H  .    i  3  <,  can  / 

jf ;  ir^  =  Tstss^/3  +  Tu;rviG*2  '3 
u;  f  =CMi/  j4j(  Jf  i  In  ) 

pk-jnf  =FNRhow '  Tf  i  Irc  ) 

rn(=:urn,,/Tf,  l  ~  i 

u  pi         .     iswp'iAniixx-.-i 
M,  ,  f  =DMM.  ,.  ./    T  P   ,    1~,\ 


\*      D   \ 


m    dc    MnnTcrrn 


wp-v-         I      i  .  vy  p  1*1   1     >      I    W  W  I   ■    ' 


Dr-,=  (Dw«l*ir7i(7i)'0/C.      -7C 

Mdvc:  =  (  'B"_0Ic£5  )_Mdv*Hsc  )/Hfg 

T  p      ADC/   (Mjjw.Mriwr   ^/Mrlwx  })     Q|       THEN 

Mrjxs/   M4wXMr4v/^    \<       C 

50T0   2593 

DMn       TD 


ll„»Mrf„»lln    /  A  v 
vv  rio.vy,,.  .-. 


T^^f.T, 


uin< 


Dr.M,       f.rnf    /L/   f 

(  I  _•  _       ^  I 

RetD"(  lJv*Do*Rhof  )/Muf 

M  .  •  m  m = I   Un*n0    \    .   L    ( 
iiuu  J  -         —    -  i  »    I 

M,,r.l»Mn../D.lr"         C 
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2639  Ff-9.81*Do*Muf*Hfg/<VvA2*Kf*Tdcf  ) 

2640  Hr,uss-.729*<Kf  A  3*9.9!  *Hfg*Rhcf A2/( Muf *Do*Tdc f  )  >A.25 
2542  PRINT 

2643  PRINT  USING  "IX,""  Uinf   Tdcf      Qfl         Ho 

2644  PRINT  USING  " ! X  ,2< DD . DD  ,2X  )  ,2( MZ . 3DE ,2X  ) " ; Uv  ,Tdc f  ,Qp ,Hc 
2£45  pojivjT 

2645  PRINT  USING  "IX,""     Nu         Retp      NuRe     Pr       F 


Hnu! 


Ucoci 


Ke; 


Hi 


Nuc 


St  CC£ f 


2549  PRINT  USING  " 1 X  ,2( MZ .3DE ,2X  )  ,2< DD.DD ,2X >  ,2( MZ .3DE  ,2X  )" ;Nusa ,Retp .Nurel ,Pr , 
Ff  ,Hnuss 

OCdA  Dmi  >M..m  /M,.,.ii 

7CC/I  a^,,c.u<  ^n.  /y,,. 

tlt.i  i-  «*  ■=  «  r  i  i..,»n .  .  u  k  ,  . ,  s  m  ,. 

i_  u    i   -r  in     ^^        i    i/lw       uiTnnuw    "    i  lulil 

•7C7C  Ar.r.««,M.....»rn... /!/... 

*_■—■_■  ■  r-  ■     ^       i  iuw'    upw    liw 

OC7C  C*nn.ar,,./(flr.,''         Q»Anr-oA         ^  ^  ~  ^  »  D  m  +    ''  1/t\ 

uu  i    i  i«upi  i  ij  ^       i  npi     _       ■  wwWfe/^rvQfa       •    l-f/ 

2g7g  ppjMT 

?CQ<7l  DDTMT      IICTMC      "IV 

LUUiy  I        I    \    A   I   M     I  OvJXIIW  I     A 

•P  M.  ,  /p- 

2£gi  ppjMT  ijsjkjc  "iv  £(NZ.3DE  2X)"'Uui  Rei  Hi  Anus  Stcc  Nucr 

2692  INPUT  "CHANGE  TCGOL  RISE7  !  ~Y  2  =  N"  Itr 

2593  IF  Itr**1  THEN  GOTO  0li\~' 

7CQ/I  nilTDIIT     «c  .   1  •  1   .  T,.,^/  *  \ 

4_   U   W  ~T  L/  W     I     I        W     I  VTI         x     .i.     ■_      .        1       I     WU    \      '■       / 

2704  jp    \~-=\    THEN 

2707  BEEP 

2710  INPUT  "OK  TO  STORE  THIS  DATA  SCT  ( • ™Y  0=N)?"  Ok s 

2713  IF  0 k a "•  1  THEN 

2725  OUTPUT  ®FilepBvcl  Bs~"  Ptran  Ft"  Enf'*)  Fr,    T1  T2  Pva^l 

"?7^g  cKin  ip 

27qc  BEEP 

2 751  INPUT  "U!IL 

9-7C/I  M»-.  ,„«  T 


TUCDC  DC  AMQTUCp  PUN   (  \   =*Y   0  =  N)r'"   I  GO 


2757  jp  j~c«=i  THEN  1777 

9770  ci  cc 

2773  IF  J-'Nrun  THEN  1777 

■?  7  o  3  c  m  n  t  p 

2933  SEEP 

2835  '  INPUT  "ENTER  PLOT  FILE  NAME"  Fdot* 

29391  AC C J  CM  @PiJe4  yn  pr,j0+* 

2542!  FOR  1=1  TO  Nrun 

?Q/ICI  CMTCD      «C,lm/l.nm      U~ 

"-w^  1 —  iv  1  »_  i  ,      £r,     xiCTturp   ,"u 


>QCC  I      MCVT      T 

.■_■'_■_!  '     "    _      .  1 

:i7ira  1        d 
-  w  *•  •        i 

n    f  .  1  -* 

t^ 1   j.  j.  c  •** 

3007      BEEP 

30101        JC       Jpf-1       TLjCM 

30  1  3         pp  J  KIT 

301£!      ppjMT     IJCJMg      "  1  0X 

""  ,10A";J ,P_file$ 
^(7i  1  a  i     cKin    tit 

-    »-  -'  UlfL/  4.   I 
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dsta  runs  were  stored  in  file  "",10A";J, 


V-V   r,.  ,  .- 


sirs  were  stored  in  plot  deto  file 
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■7. 1 


3073 
307E 

Jv  I  -J  i 

3092 

3092 
3  i  i  c 

3119 
3 1  •?  i 

3124 
■7  i  27 

~ i  33 

3E 

J9 

42 
i  4c 

/l  Q 
MCI 

54 

57 

E0 
e- 

cc 

3 1  69 
- 1  7-? 

■7  1  7c 

70 

"'•■El 
□  /i 


ACC  TCM  «C  i  1 


m 


ACC  TCM 


GN  @f 


le 


TO 


ACCTCM  «C  .  1=4.  TO   * 

1  1  _/  -/  -.  w  t  *   ci   A  x  _.  -r    1  _/ 

IF  Iso-2  THEN  CALL  Raw 

IF  Iso-3  THEN  CALL  Purg 

END 

ncc    cwpv5+ (  tc  ) 

rnM     /c  1  -4  /     t *4 

uuii       /    1     x  i_i  /         i.    1    v 

n  t  m    1/  (  o  ) 

TC      Tf  +  _.(?i      TUCM 

XI  X      I        I.  W  I     I     i  <_  I   * 

nATA      _7      CQ19^/ICC/1      -7C      p)Q£7.7  "7C  CIC      -ICO       ,7f7<CC/lF      C.  A      7  "77  Q  CC,<7M      -11Q      qp/,C77C 
nATA      /I       1CTI    1  7*77      7  (7l      q7Cfl)C7C       1      CQ     C 

U/  I    I    I     II  "T     .       I     U      I       I       I       !%_'--,_._/._Jl_jVl_ll_/,l       •    __.-,._ 

DC  An       1/  /    *    > 

t«( Tc+273. 15 )/E47.3 
c. .---«.?. 

cnD     M__ffl     Tfl     /I 

I      Wl  \       1%        _/         1    w        -+ 

Suir!"S"r!+K  (  N  )*'  1  -T  )"  (  N+  1  ) 

MCVT      M 


C-=C 


,--/  /  t»  /  1  x(/ (  c  )» (  1  _t  ) 4.1/  (  c  )*  (  1  _t  )  "9  )  )-(  1  -J  )/(  K  (  7  )*'  1  -T  )  A2+K  (  9  )  ) 


D  =  7  7  1  7l?i<?ir7i(?i»D.- 

1  *-       ^  '         A_       \-'       w        — '       _  I 

CMn    tc 

L_  I  <U  XI 

tc  t  f  +  _. i  JHEN 
Tf.j..*]  >g4.324-4cgic 

P= 1 0" ( 33 . 0E 55-4330 . 98/Tf-9 . 2S35*LST( Tf >+2.0539E-3*Tf ) 

D»D»  tf«1  79C/H      CQC 
■        i  i        -  _  -  ~»   .  l.   .  '_■ 

END    IF 

IF    I  f  t  =2    THEN 

A-9 . 394E85-30EE . 1 / ( Tc+273 .15) 
an      p=  i  ~"7    7?*i  i7i    a 

ara      CMn    tc 

_>  _/  L_  I  ■■  _  XI 

Q7  pCTIIDM      D 

qc      pjsjEND 

QQ         nCC      CMUf ->.■  T   \ 
•^  •■"- y  >     I     / 

3202   CON  /Fid/  Ift 

3205   IF  lft~0  THEN 

a7,7iq       ut~  =  i  Annr>{*n>-~>  acick*.  i  t-\  a  \ 

32 1 1       END    IF 

32:4       IF    I  f  t " 1    THEN 

3217      T f  "T*  1  .  9+"72 

*77f7l         Ufn.l     (7lCC70C7C4l_Tf»//l      C  "7Q  r?iC  7  C_  7  +  1       7C  1  ath  A  Q  C  _  A  *  T  *    , 

_/*__._.'  lliy         l,_/_J_lfw_Jl_-'l  ri\--._/_yw«_'-J--_-        _.'l._.wi->i-/Hw_.        H^ii/ 

■7777  U*„«lJ*r,»7-77C 

*_<__._;  H,y       riiy._w._u. 

-7  7C  CMTi      TC 

_■/  _—  _—  W  i -  -_-  x  I 

3229  IF    Ift~2    THEN 

3232  Tk**T  +  273.15 

3235  Hf G"1 .352E4E+E-Tk*(E.382E3E+2+Tk*.7474E2) 

a 2 3Q  ENn    IF 

-7/1  1        DCTI  idki    u;P 

_/_--rt  I\u-IWI\IN        Illy 

3244       CKjcMn 

3247      DEF    FNNi_ui(  T  ) 
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3250 
3253 

3259 

32S5 

w/  A-  W  W 

327  1 
3274 
3277 
3290 
3293 
32Q5 

3299 
3292 
3295 
3299 
330 ! 

W  W  W  l 

3310 

3313 
33  i  c 

33  1  9 
33^2 

3325 

7770 


3340 
^34"7 
3345 
3349 

~*355 
335S 

■7~C  1 

v^  w/  W      I 

3354 

3370 
3373 
3375 

77TQ 

3392 
3395 
3399 

77QI 
77Q^ 
77Q7 


rr\M      IC  1  A  I      If* 

WWII  I      \        4.1-1/  All. 

4  1  All  1U  I     I    ll_t  1 

a«247.9/(T+13'7.  15) 
Mu-2.4E-5*J0AA 

END    IF 

IF    Ift-1    THEN 

Mu-8. 962991 9E-4-T*( 1  . 1 094609E-5-T»5 . 55592 9E-9  ) 
CNin     tec 

IF    Ift-2    THEN 

TU . 1    / /  Ti?77      1  C   } 

irv  l/\l'i-l%Jml^/ 

Mm.CVDi.11       (7i  1  7D  +  TI,  +  I    1       "7  4.  /I  C  x  "S  _T  I-  *  /  7      Q0I7J7ICC  J.C-T  l  *  1       HCCIC  +  OII) 

I  i  w.        i_.-<  ii»(uiiw'in.\i»i*TTi_-w^         1    1  ■  _    .    w  u  „  w  — j  1—     ■    — ■         ir>r,r«ii_LJwiL_'w/// 

CMn       TCT 

1 —  '   ^  t —  JL  I 

RETURN   f"1" 

FNEND 

D 5 F    CMl,w3*(Tt  } 

rnivt     /ciA/    t*  + 

W  W  I    1  /I        J.    W.   /  A      I       V 

IF    Ift-0    THEN 

p=pMpv»+ ( T*    \ 

t=t++273 . 1 5 
v- 1  cpira  /T 

C1=1/(   ]+T»1   _C_4   ) 

F  2 ■  (  1  - E  X  P  (  - X ) )A2.5*EXP(X  )  / X  " . 5 

Q=      fll(7i  1  C.C  1  _      (71  (71  (71  Q  A  ?  *  C  "7-      r7ir7lf7l  /I  a  O  9  *  Y 

l;-2*P/  (  45  1  .  52^T  ) 

u  =  (  u;  i  +  2 * g » k  ) A  .5  ) / 1/ 

CKin  tc 

IF    I  ft""1    THEN 

Tf.Tt*i .9+32 

y= i 3 _ 055357- t t » ( 

u  .  u  /  ]  6  .  0 }  Q 

END    IF 

IF    I f t"2    THEN 

tl  =t-i  +273  .  1  5 

p.CKipw.^  (  J+  ) 

II  =*  1  a  7.   Q  C  *  T  p  /  D 

END    IF 

RETURN   U 

FNEND 

ncc  cwr-,.,/  t  \ 

* — '  I —  ■      I   HUfJ  W  v   I   / 

QQM      /Pljj/      Ift 

IF    lft=0    THEN 

C~u;  =  4.21  1 20959-T* 
cwn  tc 

jp    t  f  +  _  -j    JHEN 

rn..,.Q  OC(*~!TT?.C-  1  4.T*/  Q  7^(7|(7|/1  "7~C_/1x  1   V^mW^QOC-C  ,T  \ 

END    IF 

j p    t f  +  «2    THEN 

Tk~T+273 . 1 5 

Pn.,  =  /1       1QCQ*,    1       COQ1C.74.TU/7      7?C£7jQ  7C  _"7_T  I  *  ,   T      "7  "7  /t  C  _C  _T  I  *  *7     C  1  "7  4  QC  _  Q  \   S   \ 

END    IF 


TiC.Tf. 


177CI  Qf7lC_/1    -i 


7CQ7CC-7-T*/   /I       /177C1t:_Ca.?      "7  1/1-7 


A  IOC -1 *T  \ 

"T   <_    W  l_  I  '        I         / 
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•>J  T  V  \J 

7/1017: 

W  T  W  W 

7/1  01 C 

1-/  ~T  W  W 

3409 

34  1  2 

741c 

~4  1  9 

34'7 1 

3424 

342? 

3430 

-7  4  77. 

w  -r  w  w 

343B 

3439 

3442 

~744c 

3449 

345  1 

3454 

3457 

3453 

^  ~f  w  w 

3459 

3472 

3475 

-7.AIQ 

W  -T     1     W 

~49  1 

7404 

3  4  £  7 

3490 

7/97 

3495 

*499 

3502 

7C0IC 

W  W  W  W 

350S 

35 !  1 

7C  1  4 

7C  1  7 

3570 

7C_  7 

7C9C 

w  w  i_  u 

-c29 

3532 

3535 

3539 

■7.C4  1 

7C  /t  A    1 

w  w  -T  -T    : 

7C/I7 

DCTIIDM      P  n  . . ,  *  1  01 01 01 

I  \  W.     I     Wl  W  1  Up*77     I    Ut/t- 

FNEND 

pec    cmdu„v/  t  ■\ 

W  l_  I  1      I  *  IV  J   I  W  l/u    \       I       / 

rnM     /c  1  ^  /     t  *  + 

WWII  /'        AW/  1    I      I 

TC      T  P+  »0l     TUCM 

XI  X     I       V  W  I     I     I  L_  I  *l 

Rc-999.52  945+T*(  . 01 2G9-T* ( 5 . 48251 3E-3-T* 1  .234147E-5  >  ) 

CMP.   TC 

l_  .  1  W    X  1 

TC   T  f  4-  =.  1   TUCM 

4.  1         a.   1    •        I  ■   :  l_  .  1 

I  \  W     1   iUUV  I  -T  I  vJI.  '  <J  x_  .  (_  1  ■_.  U  ^  -r  □     1     L.  .  w-  —>  i  w  i_  J  I  i_    w   ' 

END  IF 

IF  Ift-2  THEN 

TK-T+273. 15-339. 15 

tif.Q     7/IQ/IQD-/I  4.TL  */  C      9  "7  DC  C --7  +  TL- *  C  Q      0  /t  /  /I  C  _  1  0,4-T  U  *  7:      (7IC7C  _  1  9  \   S 

VI  W.X,-TW-TWi_  -T     i       I     T-.  ^    \     W     .    X.     I      WW  C  I       ■       trvTr\J,X_-TfcrTW  IW'      Ir\-W.WWIL-  I    L     I      / 

R  c  "  1  /  U  f 

END  IF 

p  C  T  IJ  p  m  p  ~ 

C  M  c  m  n 

DEF  FNPru( T  ) 

D  r  w  _  C  K I |  P  r,  w  (  T  )  *  C  M  M  u  w  (  T  )  /  C  M 1/  w  (  T  \ 
p  C  T I  I O  M  D  r  y 

FNEND 

DEF  FNK'1'  T  ) 

W  W  I  '     /I    X  W  /      XI   I. 

tc  t  f * «2  TH5N 

y  -  /  y  +  9  7  3  _  i  c  )  /  2  7  3  .  1  5 

i/w=- ,  92247+ X*  >  2.8395-X*(  1  . 9007-X* (  . 52577-. 07344*X  )  )  > 

END    IF 

tc    if*i]     THEN 

Ku,,  =  c    20952"r9E-2-T*(2    ~",M~'95E-4  +  T'*•',    38095245-8) 
c  m  n    t  c 

IF    I  ft  =  2    THEN 

ti =T+273 . 1 5 

Koj  =  4.  1859E-4*(519.442+.320920*Tr;  ) 

ENn    IF 

RETURN   Kw 

FNEND 

ncc    cwT,rk( y  \ 

WL_I  I        I    *      I      W   f    I    I     I     \      *\       / 

P-EXP( X  ) 

n  =  cvD/ _v  \ 

t.^u./  D_n  n ii  Din  \ 

DCTIIDM      T  =  ~U 

I    V   l_       I       W  I    V  1    S  I      W    I     I    I     I 

CM  CMP) 

nC_C     CMTw.  v  (  IJ  ) 

WW1!    I         -     WW   -  W    v      I      /  , 

enp     T  =  1     tq    7 

T  "  T  +  C  (  I  )  *  lJ  A  I 
NEXT    I 

T.T  +  /1      7770CC-7iTx(7      CQTQ74C-7-T.D      (7I"77  q  *?  "7C  _C  ) 

I  |l    -T.l^/WWWL-  W'l^VI.WWWWW-TW.  W  I^W.WIIWX-IW  Wl 

DCTIIDM      T 
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"3  c  ->  2  p  <-  -  p  w  p  v  5  t  '  T  "  ) 

3525  jp  ftgc/ / P-Pc )/P )> . 0001  THEN 

3E2G  IF  Pc'P  THEN  Ti=T3 

353;  jp  Pc>p  THEN  Tu^Ts 


3550 
3553 

355E 
3559 
3552 
3555 

7CCD 

357  1 
3574 
3577 

^J   wl  w  t/ 

35C3 

7COC 

w/  U  W  Ul 

3599 
3552 

7CQC 

W  *J  *J  w> 

35  i  g 
35  i  3 

7C  1  C 

W  W  I  t-J 

-%C   1  Q 

a5">2 
3525 

553  i 

3c~r'4 
3537 

3540 

3543 

39  5  3 
755- 

7C77 
7.CD7 

w'   W   W   ~> 

■3"7m^ 
37 1  3 

7"i~ 
"a*7"7  "7 

3734 

3735 
3735 

~7",:7 
3747 
3749 
3749 
3753 

7TC/I 

3755 
3755 

3757 

77C7 


IE0317E-E ) ) 


FNEND 

QCp  pMLJf  (  J  ) 

COM  /Fid/  If t 
IF  Ift-0  THEN 
Hf-T*(4.203849-T*<5.88132E-4-T*4.551603 

END  IF 

TC   TP*»I   TUCM 

XI       Alt      I      I  I  ll_  l\ 

Tf "J* 1 . 9+32 

Hf-8. 207957 1+Tf*(  . 1 94E7957  +  Tf *  1 .321429EE-4 ) 

Uf.Uf *.0      79C 

III  ill  __        .        -_      —      U 

END    IF 

IF    Ift~2    THEN 

Uf-7C(7l    |   Tn  DC   UCDTpTCn 

III       ^  w'  L  I   W     UL      */  I—  I    A  I    J.  l_  t^ 

END  IF 

pCTIJDM  U  f  *  1  £  0  0 

FNEND 

net:  :wTu.  „  /  p  \ 

U/t-l      I   I  *  I   V  O  ^  \  I    / 

Tu" 1 90 
Tl  - 1  0 

Ts~(Tu  +  TI  )  *  .  5 

Dr.CMDw. t  (  T' 


GOTO  35 ! 9 
END  IF 

RETURN  Ts 

FNEND 

n  r  r  p  m  t  v  s  -y-  5  5  (  lJ  ) 

COM  'CcS1^/  TF5(5) 

T  =  T55( 0  ) 

FOR  1=1  TO  5 

T*T  +  T55(  I  )*U"  I 

NEXT  I 

RETURN  T 

FNEND 

DEF   FNTv3v5E(v' ) 
rnivt     /r~cc/    tcc/c\ 

w     W     '       I  -      -—      —I    LJ  >        _■!_'-        w  ' 

T.TCCf  tfl   \ 

■  i    _»  U    v    \J     i 

pnp    t=i     yn    c 

7  =  74.755  (   T    )*U  "  T 

NEXT    I 

RETURN    T 

FNEND 

npp  FNTvsv57(  'J  ) 

COM  /Cc57/  T57(5) 

T~T57( 0  ) 

pnp  j«]  7Q  5 

T~T+T57( I )*U" I 

MtTVT   T 
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37E8 
3773 

-T7A 

*J    i  I  -r 

3775 
->.  -77c 

■3777 

A-7D-7 

3797 
■?ora-7 

-  -  *~    1 

7Q  1  "7 

W  W  1   I 

3827 

7Q77 

3847 

7QCT 

W  W  w>  I 

38E7 
3877 

7Q7D 

3879 

^OQr7l 

_  _  _  1, 

7QO] 

3883 
'  984 
3885 

7QOC 

w  \^  w  w 

3887 
3888 
398° 
3890 
389 1 

ICG'' 

WW    Jl. 

3993 

aQQ/1 

3895 

7oac 

_>  ^  -  - 

■7007 

3998 

3899 

~Qf7ll7l 

w'   W  IL/  *^ 

3902 
3903 
3904 
3906 
3907 
3909 

3917 

7Q?7 
7D77 


RETURN  T 
p  m  rwn 

nCC  pMTw«vCQ(  U  ) 

COM  /Cc58/  T58( 5  ) 

T.TCQf  a>  \ 
pnp  t_i  yn  p. 

T.TxTCQ(  T  N  +  l  J'-  T 

NEXT  I 

RETURN  T 

FNEND 

DEF  FNCcsh( X  ) 

p  ■  c  v  p  (  V  ) 

Q-EXP(-X  ) 
Cosh"  c* 1 P+n  ) 

D  C  T  I J  D  |\|   P  ~  a  J-i 

p  m  d  m  n 
SUE  R  3  u 

POM   /CU/   Tf* 

\-<  W  I  1    /I   1  U'      0.1   V 

DIM  Enf  (  13) 
BEEP 

T  MDI  IT 

J.  I  <t  I   W  I 

ASSIGN  @Fi 
BEEP 

T  MDI  IT 

4.  I  "4  I   W  I 

INPUT 


cmtcd  pj_[_p  NAME 


dmtdd  pppccijpp  CONDITION  (  0  =  U 
rwTcp  NUMBER  OF  RUNS"  Nrjn 


1  =A   \ 


Tr^ 


D  D  T  M  T  D  D      TC      701 
DD  TMT 

1      <   \   A  I  V    I 


PRINT 

DDTMT 


IC  TMC 


Tin 


TOLil 


TC      Tr-,-«=f7l      TUDM 

XI  j.  H   U  <J  f    I    ,1-11 

ppTMj  US  INC  "  1  0X  ""Pre  5  5  Li  re  Cor,  di  t  ion  :  Uscui"n 

ELSE 

DDTMT     IICTKIC      "  1  0IV      "  »  D  ^.»  =  =  .  ,  r-  «>     n^^rH  +  ion'      A  +  .n~  = -*^«».-< 

I        l\iH     I  -      -       x         (    _  V     ■  .  I        I         ____"■        '_  WUI    IU1      i.     1UI    I     •  ll     VI    IW   JjJIIUI 

END    IF 

DM  TDD     RCiU.Tfn      T„„ 

DM-j-pp   @File*Inujt    Fd    Fu;   Fh 
PRINT 

DDTMT     IIC  TMC      "  1  |7|V       »  »  n  .  +  » 

DQp      To]       TQ     Mrijri 

DMTDD      flC.lp.Dw^l       D_~,r,      D  4-  „  _,  . 

DDTMT      D.,^1       D«~,r,      D  +  -  =  -      C  +  ~      r~f  C   «   1      D~>      T1       T7      D 

NEXT    I 

ACCTCM     «D 

SUBEND 

CIID      D..-~ 
^wE     I    Li i    ^ 

DDDD 

DLLi 

T  MDI  IT      "DMTDD      DTI    D     MAMD      TO     DD      nri    DTDO"      D,    1«* 

xillwi  ^    •   '     ■     —  I        il.L         I'll    II    It—  I     W  L>  I—  L.    i-   I u_    u_  ,l        j.      >.     _    ■*" 

PURGE  F: I eS 
GOTO  3908 

CI  IDDKin 

JUUI-llU 


,r-n 


T  1 


T"7 
>n  1 


D, 


.e    TO    * 
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